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Chapter  1 
INTRODUCTION 


An  important  topic  in  electromagnetic  field  analysis  is  the 
analysis  of  waveguide  discontinuities.  This  is  because  these  dis¬ 
continuities  are  basic  components  in  many  microwave  devices*  so  that 
an  exact  knowledge  of  their  characteristics  is  essential.  Also,  the 
analysis  of  waveguide  discontinuities,  from  a  theoretical  point  of 
view,  is  an  interesting  problem.  The  theory  developed  has  found 
immediate  application  to  some  types  of  discontinuities,  and  has  in¬ 
deed  deepened  the  understanding  of  electromagnetics  in  general. 

Despite  the  extensive  consideration  of  the  subject,  however, 
the  analysis  has  been  confined  to  a  small  set  of  discontinuities. 
Specifically,  discontinuities  that  are  coincident  with  or  have  regu¬ 
lar  cross  sections  symmetric  about  one  waveguide  cross  section  and 
are  uniform  along  one  of  the  waveguide  axes  have  been  the  most  fre¬ 
quently  treated  ones.  Furthermore,  in  a  typical  situation,  only  a 
limited  number  of  discontinuities,  usually  one  or  two,  of  the  same 
geometry  would  have  been  assumed.  General  systems  of  discontinuities, 
it  appears,  have  never  been  considered  before. 

In  analyzing  waveguide  discontinuities,  the  complete  field 
solution  is  of  very  little  interest.  Rather,  the  effect  of  the  dis¬ 
continuities  on  the  incident  modes  is  what  must  be  described  as 
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accurately  as  possible.  From  an  engineering  perspective,  descrip¬ 
tions  employing  networks  of  lumped  elements  are  preferred.  Any  such 
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network  must,  for  a  structure  that  is  both  lossless  apd  reciprocal, 
obey  the  two  basic  network  laws:  the  conservation  of  complex  power  law 
and  the  reciprocity  law. 

One  pf  the  earliest  methods  used  for  determining  the  elements 
of  the  network  representation  of  waveguide  discontinuities  has  been 
the  "Variational  Method"  [1].  In  this  method,  the  network  elements 
are  so  expressed  that  they  are  stationary  with  respect  to  arbitrary 
small  variations  of  some  field  related  quantity  about  its  true  value. 
That  is,  the  network  elements  would  be  accurate  to  an  order  one  higher 
than  that  of  the  trial  quantity,  A  judicious  choice  of  the  trial 
quantity  can  then  lead  to  remarkably  accurate  results. 

Although  ingenious  and  powerful,  the  application  of  the 
Variational  Method  has  primarily  been  restricted  to  single  discon¬ 
tinuities  of  simple  geometries  and  moderate  sizes  that  are  symmetric 
about  a  waveguide  cross  section.  This  is  because  finding  out  the 
appropriate  trial  quantity,  except  perhaps  for  some  simple  geome¬ 
tries,  requires  a  great  deal  of  insight  into  the  problem,  and,  at 
times,  solving  another  problem.  Furthermore,  although  the  results 
obtained  can  be  improved  in  a  systematic  manner  using  standard  pro¬ 
cedure  [2,  Section  7-6],  the  process  is  quite. laborious .  An  exten¬ 
sive  collection  of  theoretical  and  numerical  results  for  a  large 
variety  of  waveguide  discontinuities  obtained  using  this  method  can 
be  found  in  the  Waveguide  Handbook  [3], 


The  Variational  Method  has  gradually  given  way  to  the  more 
general  viewpoint  of  the  "Method  of  Moments11  [4]  •  Because  of  its 
variational  character  [5],  and  combined  with  the  very  fast  develop¬ 
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ment  of  computer  systems  and  software  techniques,  the  Method  of 
Moments  has  become  a  basic  tool  in  the  study  of  waveguide  discon¬ 
tinuities,  as  well  as  of  many  other  areas  of  electromagnetics.  Here, 
the  network  representation  is  determined  in  terms  of  some  field  re¬ 
lated  quantity,  for  which  an  operator,  usually  integral,  equation 
is  to  be  solved.  The  solution  proceeds  by  expanding  the  unknown  field 
quantity  as  a  linear  combination  of  some  appropriate  functions.  En¬ 
forcing  the  governing  equation  in  some  way  then  leads  to  a  matrix 
equation  for  the  coefficients  of  the  expansion,  which  can  be  solved 
using  a  matrix  solution  routine  [6]. 

Moment  solutions  are  analytically  simple*  iti  the  sense  that 
no  excessive  manipulations  are  needed  in  order  to  minimize  the  compu¬ 
tational  phase  of  the  solution.  Furthermore,  the  amount  of  effort 
expended  in  solving  for  many  discontinuities  is  the  same  as  that  is 
put  solving  for  only  one.  Since  computer  codes  are  always  written, 
the  automatic  improvement  of  the  solutions  is  very  easily  done.  The 
Method  of  Moments,  therefore,  is  well  suited  for  solving  systems  of 
waveguide  discontinuities,  provided  a  suitable  analysis  can  be  given. 
Recently,  the  application  of  the  Method  of  Moments  to  waveguide 
discontinuities  has  been  an  area  of  active  research.  Specifically, 
single  and  triple  inductive  posts  of  circular  cross  section  have 
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been  considered  in  [7],  [8],  where  many  of  their  characteristics 
have  been  discussed. 

As  is  pointed  out  earlier ,  general  systems  of  waveguide  dis¬ 
continuities  have  yet  to  be  considered.  It  is  the  purpose  of  this 
dissertation  to  consider  three  such  systems. 

The  first  system  is  that  of  multiple  apertures  of  arbitrary 
shape  in  the  transverse  plane  between  two  cylindrical  waveguides. 

The  second  system  consists  of  metallic  obstacles  in  a  rectangular 
waveguide  that  are  uniform  along  the  narrow  side  of  the  waveguide, 
but  are  otherwise  of  arbitrary  shape  and  thickness,  i.e.,  a  system 
of  inductive  posts.  The  third  system  consists  of  metallic  obstacles 
in  a  rectangular  waveguide  that  are  uniform  along  the  broad  side  of 
the  waveguide,  but  are  otherwise  of  arbitrary  shape  and  thickness, 
i.e.,  a  system  of  capacitive  posts.  Common  between  the  first  and 
second  systems  are  the  inductive  windows  and  strips  in  a  rectangular 
waveguide,  and  between  the  first  and  third  systems  are  the  capacitive 
windows  and  strips  in  a  rectangular  waveguide. 

The  system  of  multiple  apertures  is  considered  in  Chapter  2, 
and  is  depicted  in  Figure  1  there.  Assuming  a  multi-mode  exciting 
field,  a  modal  expansion  is  used  to  express  the  field  in  the  two 
waveguides.  A  field  equivalence  theorem  and  Galerkin  procedure  are 
then  utilized  to  obtain  the  generalized  network  representation  of  the 
apertures.  This  representation  is  shown  to  obey  the  two  basic  net¬ 
work  laws:  the  conservation  of  complex  power  law  and  the  reciprocity 
law.  The  scattering  matrix  is  then  deduced  from  the  generalized 
network  representation,  and  its  properties  are  examined.  The  analysis 
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is  subsequently  specialized  to  the  problems  of  inductive  and  capaci¬ 
tive  windows  in  a  rectangular  waveguide,  where  the  dominant  mode  is 
the  only  incident  wave.  The  impedance  matrices  of  the  windows  are 
also  obtained,  and  readily  realized  by  networks  of  shunt  reactive 
elements . 

The  system  of  inductive  posts,  which  can  be  seen  depicted  in 
Figure  1  of  Chapter  3,  is  then  considered.  The  exciting  field  is 
taken  to  be  the  dominant  waveguide  mode.  A  complete  field  analysis 
of  the  problem  is  given.  The  analysis  is  quite  general,  and  results 
in  an  integral  equation  for  the  currents  induced  on  the  posts.  Later, 
this  equation  is  solved  through  a  Galerkin  procedure.  The  scattering 
and  impedance  matrices  describing  the  effect  of  the  posts  on  the 
dominant  waveguide  mode  are  then  obtained  and  examined  in  detail .  The 
impedance  matrix  is  subsequently  realized  by  a  two-port  T-network  of 
reactive  elements.  The  computed  reactances  of  some  post  configurations 
are  also  reported. 

The  system  of  capacitive  posts  is  then  considered  in  Chapter  4, 
where  it  is  depicted  in  Figure  1.  Since  this  system  is  dual  to  that 
of  the  inductive  posts,  the  analysis  is  drawn  on  similar  lines.  Spe¬ 
cifically,  a  complete  field  analysis  leads  to  an  integral  equation 
for  the  currents  induced  on  the  posts  due  to  an  incident  dominant  mode 
wave,  which  is  later  solved  through  a  Galerkin  procedure.  The  scat¬ 
tering  and  impedance  matrix  representations  of  the  posts  are  then 
extracted  from  the  analysis,  and  their  properties  are  discussed.  The 
impedance  matrix  is  subsequently  realized  by  a  two— port  T— network  of 
reactive  elements.  The  computed  reactances  of  some  selected  post 


configurations  are  given. 
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A  word  about  the  organization  of  the  dissertation  is  in  order. 
Each  system  of  discontinuities  is  considered  completely  independent 
of  the  others.  This  is  believed  to  be  the  best  approach  to  the  sub¬ 
ject,  despite  the  repetition  of  some  of  the  definitions  and  proofs. 

A  discussion  at  the  end  of  each  chapter  points  out  the  important 
results  of  the  analysis  and  closely  related  works.  Final  remarks  are 


given  in  Chapter  5. 
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Chapter  2 

MULTIPLE  APERTURES  IN  THE  TRANSVERSE  PLANE 
BETWEEN  TWO  CYLINDRICAL  WAVEGUIDES 


Consider  a  system  of  apertures  S  ,  S  ,  ...,S^  of  arbitrary 
shape  located  in  the  transverse  plane  between  two  uniform  cylindrical 
waveguides  A  and  B  extending  along  the  z-axis.  The  mediums  filling 
waveguides  A  and  B  are  assumed  linear,  homogeneous,  isotropic,  and 
dissipation  free,  and  are  therefore  characterized  by  the  real  scalar 
permittivities  and  and  the  real  scalar  permeabilities  y^  and  y^> 
respectively*  Figure  1  shows  the  problem  at  hand. 


1 .  Basic  Formulation 

Let  a  multi-mode  field  be  incident  in  waveguide  A.  Part  of 
the  incident  field  is  then  reflected  into  waveguide  A,  while  the  rest 
of  it  is  transmitted  into  waveguide  B. 

The  total  z-transverse  field  in  both  waveguides  can  be  ex¬ 
pressed  in  modal  form  as  [2,  Section  8-2] 


E  = 


_ 


Yai2 


-y  .  z 

Yc.e  31  e  .  +  Y  a.  e  e  z  <  0 

h  i  — ai  h  1  — ai 
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P'YbiZ 
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(1) 
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—  — bi 


z  >  0 . 
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Figure  1.  Waveguides  A  and  B  opening  into  each 
through  S1,S2,...,SP. 


other 
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All  the  modes  TE  and  TM  to  z  are  included  in  the  summation.  In  (1)5 
ci»  a^,  and  b^,  are  the  amplitudes  of  the  ith  incident,  reflected, 
and  transmitted  modes,  respectively.  Ya^  and  r|  ^  are,  respectively, 
the  modal  propagation  constant  and  characteristic  admittance  of  the 
ith  mode  in  waveguide  A: 


ai 


X  <  X  . 
a  ai 

(2) 


X 

a 


>  X 


ai 


Yai 

jwu 


Yai 


for  TE  to  z  modes 


for  TM  to  z  modes. 


(3) 


Here,  K  is  the  wave  number  of  the  waveguide  medium,  k  .  is  the  ith 
a  ax 

mode  cutoff  wave  number,  and  X  and  X  .  are  the  corresponding  wave 

a  a  x 

lengths.  The  parameters  for  waveguide  B,  and  are  similarly 

defined.  Finally,  the  modal  vectors,  e^  in  waveguide  A  and  in 
waveguide  B,  are  assumed  real  and  so  normalized  that 


e  .  .  e  ,  ds 

J  *qi  -qj 
Q 


i  -  1 
i  +  j  • 


(4) 


The  integration  in  (4)  is  taken  over  the  cross  section  of  waveguide 
Q  e  {A,B} ,  and  q  £  {a,b}  is  such  that  (Q,q)  ~  (A, a)  or  (B,b). 

For  a  complete  field  solution,  the  amplitudes  a^  and  b^  of 
the  various  modes  in  the  two  waveguides  are  to  be  determined.  This 
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can  be  accomplished  with  the  help  of  a  field  equivalence  theorem. 

Let  the  exciting  field  be  incident  in  waveguide  A  while  Sm  are 
covered  by  perfect  conductors.  This  field,  sometimes  referred  to  as 
the  generator  field,  is  denoted  (Er,  H  ).  The  field  equivalence 
theorem  [2,  Section  3-5]  states  that  the  field  in  waveguide  A  is  iden¬ 
tical  with  (.Eg,  H®)  plus  the  field  produced  by  the  magnetic  current 

sheet  M  =  M™  where 
m=l 


l  X 


on  S 


m 


(5) 


while  each  S™  is  covered  by  a  perfect  conductor.  The  field  in  waveguide 
B  is  then  identical  with  the  field  produced  by  the  magnetic  current  -M 
while  each  S™  is  covered  by  a  perfect  conductor.  Figure  2  shows  the 
equivalent  situations . 

The  z-transverse  field  produced  by  M  in  waveguide  A,  denoted 

(E  (M)  ,  H  (M) ) ,  and  that  produced  in  waveguide  B  by  -M,  denoted 
—a  —  —a  — 

(E,(-M)  ,  EL  (-M) ) ,  have  the  same  form  as  (1),  except  that  there  is  no 
_b  —  — b  — 

exciting  field.  Thus,  the  total  z-transverse  field  in  both  waveguides 
is  given  by 


-y  .z 

'ai 


y  .z 
!ai 


y  .  z 
’ax 


E^  +  E  (M)  =  J1  c.e  e.-Yc.e  e.+Vd.e  e. 
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(a) 


Figure  2.  (a)  The  equivalent  situation  for  waveguide  A. 

M  exists  only  on  S  ♦ 

(b)  The  equivalent  situation  for  waveguide  B. 
_M  exists  only  on  S  , 
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In  (6),  c^,  d^,  and  b^  are,  respectively,  the  amplitudes  of  the  ith 
incident  mode,  the  ith  mode  produced  in  waveguide  A  by  M,  and  the  ith 
mode  produced  in  waveguide  B  by  -M.  It  then  follows  from  (5)  and 


(6)  that 


M  =  < 


V  d .  z  x  e  .  on  A 

l  ±  —  -ai 


z=0 


y  b .  zxe 

l  x  —  — bi 


on  B 


z=0 


(7) 


The  placement  of  magnetic  current  sheets  Mm  in  waveguide  A  and  -Mm 
in  waveguide  B  over  each  Sm  therefore  ensures  the  continuity  of  12 
across  the  apertures.  The  continuity  of  _Ht>  however,  requires  that 


2  ? 


c  .n  .  zxe  . 
l  ai  —  — ai 


.zxe. 
ai  —  — ai 


on 


m 


+  I  b.n.  ,  zxe,. 
h  i  bi  -  — bx 
x 

1  <  m  <  p 


(8) 


which  is,  together  with  (7),  the  equation  determining  the  amplitudes 

d .  and  b . . 
l  1 
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2 .  The  Generalized  Network  Representation 


An  exact  solution  of  (7)  and  (8)  for  the  amplitudes  d^  and  b^, 
and  consequently  the  complete  field  solution,  can  rarely  be  obtained. 
However,  only  a  representation  of  the  apertures  that  describes  their 
effect  on  the  modes  of  the  waveguides,  not  the  amplitudes  of  these 
modes,  is  usually  all  that  is  needed.  In  this  section,  a  representa¬ 
tion  in  terms  of  two  generalized  networks  is  derived. 

Since  the  set  {e^!  1  -  1,2,...}  is  complete  [9,  Section  5-6], 
a  finite  subset  of  the  lower  order  modes  can  be  used  to  approximate 

the  field  in  waveguide  Q.  Thus,  (7)  and  (8)  become 

£ 

a 


l 

i=l 


b 

l 

i=l 


d,  z  x  e  . 
l  —  — ai 


b. 

l 


zxe.  , 
—  — bi 


on 


A 


z=0 


(9) 


on 


B 


z=0 


i 


2  y  c.T]  .  zxe  ,  =  y  d.n  ,  zxe  ,  +  J  b.rt.zxe 

i  ai  —  — ai  ^  l  ai  —  -ai  .  l  bi  —  — bi 

i=l  i=l  1=1 


m 

on  S  ,  1  <  m  <_  p  (10) 


where  £  and  JL  are,  respectively,  the  number  of  modes  used  in  the 
a  b 

modal  expansion  of  the  fields  in  waveguides  A  and  B.  Here,  d^  and  b^ 

are  no  longer  the  actual  mode  amplitudes.  Only  in  the  limit,  as  £^ 

and  JL  go  to  infinity,  they  become  so.  However,  for  sufficiently 
b 

large  £  and  £  ,  it  can  be  assumed  that  d.  and  b.  in  (9)  and  (10) 
a  b  11 

represent  the  actual  mode  amplitudes,  while  equality  still  holds 


there . 
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Let  (M;  |l  <.  j  p  }  be  a  set  of  real  vectors,  and  put 


m 


p  m  m 

„m  r  V.  M. 
M  =  L  J  ~3 
j=l 


(ID 


where  V™  are  complex  coefficients  to  be  determined.  Substituting 

(11)  into  (9),  it  becomes 

£ 

a 

V  d .  z  x  e  .  on  A 

. l  —  — ai 

i=l 


m 


l  Vm  M .  = 

j=l  3  ~3 


l 


b 

y  b .  z  x  e,  .  on  B 
. l  —  — bi 

i=1 


z-0 


z=0 


(12) 


Put 


w 


<lk  . 


M  •  z  x  e  _  ds  . 
-  -  -qk 


z=0 


(13) 


It  then  follows  from  (4)  and  (12)  that 


»ak-  l  f  vjj 

ak  m=l  j=l  3  J 


Mm  •  z  x  e  ,  ds  =  d, 
— j  —  — ak  k 


m 


%  r 

=  l  b. 

1=1  1 


zxe,  •  z  x  e  ,  ds,  1  <  k  <  £ 
—  — bi  —  — ak  —  —  a 


(14) 


aRb 


m 


v  V  m 

Wbk  ~  t  .  -i  Vj  J_  — 3  ~  "  ~bk 

m^l  j=l  im 


M*!*  •  z  x  e.,  ds 


=  I  d. 

1=1  1  i 

11  AflB 


zxe  .  •  z  x  e, ,  ds  =  b,  ,  1  <  k  <  JL  • 

—  — ai  —  — bk  k  —  —  b 


(15) 


Finally,  scalarly  multiplying  (10)  by  and  integrating  over  Sm, 


it  becomes 
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2  E  i 

1=1  1  ai  J 


M™  •  z  x  e  . ds  =  I  d . n 
m  -3  -  ~ai  i  ai 


m 


M.  •  z  x  e  .  as 
-j  -  -ai 


£k 

b  r 

+  ^  binbi  MT  *  1  x  £bids»  1  £  m  <  p,  1  <  j  <  p”  . 
i=!  J,  m 

(16) 


In  matrix  form,  (14),  (15),  and  (16)  become 


-V  —b  *► 

w  =  H  V  =  d  =  H  b 
a  a 


w,  =  H,V  =  Hd  =  b 
b  b 


T  -*  T  ->■  T  -* 

2  HI  c  =  HI  d  +  H^Y  ,  b  . 

a  oa  a  oa  b  ob 


(17) 

(18) 

(1?) 


Here,  T  denotes  matrix  transpose,  w  is  the  vector 


w  =  [w  . ]  n 
q  q3  £q  *  1 


(20) 


V  is  the  p  segment  vector  whose  mth  segment  is  the  vector 


vm  =  [vm] 

J  pmxl 


(21) 


c,  d,  and  b  are  the  vectors 


c  [cjh  xi 

J  a 


*-IdAxi 


b  = 


(22) 

(23) 

(24) 


H  is  the  p  block  row  matrix  whose  mth  block  is  the  matrix 

q 
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Hm  = 


[H®  J 


£  xp11 

q 


=  [ 


’  m 


z  x  e  .  •  M,  ds] 

-  "J  £xpra 

q 


(25) 


Y  is  the  diagonal  matrix 
oq 


^oq  ^oqii^£  *£  ^qi^£  x£ 

q  q  q  q 


(26) 


and  H  is  the  matrix 


H  -  ("ij'vi 
J  b  a 


ARB 


£  *  £ bi  •  i  *  £ajdslt  XU  • 

b  a 


(27) 


Since  H  ,  H,  ,  and  H  depend  only  on  the  sets  of  waveguide  modes  used 
a  b 


I  p  I  m 

in  the  modal  expansion,  the  set  SV  (m‘“}  of  magnetic  currents, 

m=l  3=1  — j 

and  the  shape  of  the  apertures,  it  readily  follows  from  (17)  and 
(18)  that 


H,  =  HH 

T  b  s 


T 

HH  =  U 


HH  =  TJ  . 


(28) 

(29) 

(30) 


U  is  the  identity  matrix  of  order  £g  in  (29),  and  of  order  £fe  in 
(30).  Using  (17)  and  (18),  (19)  becomes 


(Y  +  Y,  )V  =  I 
a  b 


(31) 


where 


Y  =  H  Y  H 
q  q  oq  q 


(32) 


T  -* 

I  =  2  H  Y  c 
a  oa 


(33) 
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By  (31) ,  the  generalized  network  representation  of  the  aper¬ 
tures  is  finally  obtained.  This  representation  consists  of  two  net¬ 
works  Y  and  Y,  in  parallel  with  the  current  source  I,  as  (31) 
a  b 

readily  indicates.  Figure  3-  depicts  such  a  representation.  In  the 
following  two  sections,  the  generalized  network  representation  is 
shown  to  obey  the  two  basic  network  laws:  the  conservation  of  com¬ 
plex  power  law  and  the  reciprocity  law. 


3 .  The  conservation  of  Complex  Power  Law 

The  complex  power  transmitted  through  the  apertures  into 
waveguide  B  is  basically 


tr 


E  x  H  •  z  ds 


(34) 


B 


z=0 


where  *  denotes  complex  conjugate, 
and  using  (13),  (34)  becomes 


Substituting  from  (5)  and  (6), 


L 


p  .  i  b* 
tr  i-i  1 


*  * 


nbi  wbi 


(35) 


In  matrix  form,  (35)  becomes 


-*H  *  -»• 

P„  =  b  Y  ,  w, 
tr  ob  b 


(36) 


In  (36),  H  denotes  matrix  Hermitian.  It  then  follows  from  (18) 
and  (32)  that 


tr 


=  b  V  b 

ob 


VH  H**Y*  h  V  =  vV  V 
b  ob  b  b 


(37) 


Thus,  the  complex  power  transmitted  into  waveguide  B  is  equal  to  the 


-L. 

1  H 

(t> 

v,1 

SB 

Y0 

«P 

Yb 

+  i 

- 

<!> 

+  D 

.  Pp 

Figure  3.  The  generalized  network  representation  of  the 


apertures . 
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complex  power  at  the  terminals  of  network  Y^. 

The  total  complex  power  entering  the  apertures  from  waveguide 
A  is  given  by 


in 


E  x  H  •  z  ds 


z=0 


Substituting  from  (5)  and  (6),  and  using  (13),  (38)  becomes 


(38) 


£ 

r  *  *  * 

P.  =  >  (2c.  -  d  )n  .w  . 

m  l  i  ai  ai 

i=l 


(39) 


In  matrix  form 


P.  =  2c  Y  w  -  d  Y  w  . 
in  oa  a  oa  a 


(40) 


It  then  follows  from  (17),  (32),  and  (33)  that 


P.  =  2c 
in  oa 


Y*  d  -  d^Y*  d  =  2cHY*  H  V  -  vV  Y*  H  V 


oa 


oa  a 


a  oa  a 


-  n  -  v  . 

a 


(41) 


Thus,  the  total  complex  power  entering  the  apertures  from  wave- 

guide  A  is  equal  to  the  complex  power  supplied  by  the  source  I  minus 

that  at  the  terminals  of  network  Y  . 

a 

The  conservation  of  power  law  states  that  the  total  complex 
power  in  any  network  must  be  zero,  or 


P  -  P  =  0  . 
tr  in 


(42) 


This  is  guaranteed  because  of  (31).  The  conservation  of  complex 
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power  law  holds  if  and  only  if  the  complex  power  flow  across  the  aper¬ 
tures  is  continuous. 


4 .  The  Reciprocity  Law 

Lei  a  multi-mode  field  be  incident  in  waveguide  B.  The  total 
z-transverse  field  in  both  waveguides  is  then 


Ef  = 
—t 


Y  -z 

r  ji  ai 
)  d  e  e  . 
r  i  — ai 

l 


z  <  0 


Y,  .  2  “Yu  * 2 

V  c 1  e  ^  e  +  y  a '  e  1  e...  z  >  Q 

h  i  — bi  “  i 

i  i 


-bi 


(43) 


HT  =  < 


Y  -2 

-  y  d!n  .  e  ai  z  x  e  , 
V  i  ai  —  — ai 

x 


Ybiz 


z  <  0 


'Ybiz 


z  x  e. 


-  V  c  n,  .  e  bi  z  x  e,  .  +  I  a'n,  .  e  _ 

i  'bi  ~  — bi  r  i  bi  -  ~bi 

l  i 

z  >  0 

where  c!,  a!,  and  d!  are,  respectively,  the  amplitudes  of  the  ith 
ill 

incident,  reflected,  and  transmitted  modes. 

The  generator  field  (E'S,  H'8)  is  the  field  that  would  exist 
if  the  exciting  field  was  incident  in  waveguide  B  while  all  S™  were 
covered  by  perfect  conductors.  By  the  field  equivalence  theorem,  the 

o-  cr 

z— transverse  field  in  waveguide  B  is  identical  with  (E/  9  )  plus 

the  z-transverse  field  (E^(-M*),  H£(-M'))  produced  by  the  magnetic 

current  sheet  -M!  =  -  M™  where 

m=l 


M'™  =  z  x 


on  S 


m 


(44) 
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while  each  Sm  is  covered  by  a  perfect  conductor.  The  z-transverse 
field  in  waveguide  A  is  then  identical  with  the  z-transverse  field 
(_E'(M')}  H'(M'))  produced  by  the  magnetic  current  M'  while  each  S™ 
is  covered  by  a  perfect  conductor. 

The  total  z-transverse  field  in  both  waveguides  is  then 


K 


Y  .  z 

E  *  (M ' )  =  £  d !  e  e, 
-a  —  t"  1  ~ai 


z  <  0 


IJg+E^(-M')  =  l  cje  bl  £bl-  l  cje  bl  ebl+  l  b!e  bl  e 
i  i  i 


-bi 


-  < 


Y  .2 

=  -  l  d^nal  e  z  x  e 
i 


Ybiz 


ai 


z  >  0 


z  <  0 


(45) 


~YbiZ 


IbS+H^(-M')  =-  I  c^nbi  e  z  x  £bi  -  l  c|nbie  z  x  e 


-bi 


-Yhiz 

+  y  b!n,  .  e  z  X  e 

r  i  bx  — 


bi 


z  >  0  . 


In  (45),  c! ,  d^,  and  b^  are  the  amplitudes  of  the  ith  incident  mode, 
the  ith  mode  produced  by  Mf  in  waveguide  A,  and  the  ith  mode  produced 
by  -M?  in  waveguide  B,  respectively.  The  continuity  of  across  the 
apertures  implies  that 


M 1  = 


V  d1  z  x  e  .  on  A 

Li  —  —ai 


h  i 
i 


Y  b!  z  x  p  on  B 

l  x  —  — bx 


^  1 


z=0 


z=0 


(46) 


t 


whereas  for  that  of  H1 
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2  5!  Q !  n.  .  z  x  °  y  d'r)  .  z  x  e  ,  +  y  b'n,  .  z  *  e, 

J  i  bi  -  — bl  J  i  'ai  -  -ai  J  i  'hi  *  -M 


„m 


on  S  ,  1  <  tn 


(47) 


must  be  satisfied. 

Let  a  finitp  subset  of  the  lower  modes  in  waveguide  Q  be  used 
tq  approximate  the  field  there.  The  cardinality  of  this  set  is  set 
equal  to  A  .  Put 
m 

M,tn  =  \  V!m  M®  (48) 

j»l  J  J 


wqi  "  j  M'  •  z  x  e^  ds  (49) 

^  |  z*=0 

where  y^m  are  complex  coefficients  to  be  determined.  This  situation 
is  analogous  to  that  in  Section  2,  and  can  similarly  be  treated .  Thus, 
in  analogy  with  (17),  (18),  and  (19),  one  ohtains 


w'  *=  H  V’  -  d*  »  HTb'  (50) 

a  a 

v'  **  H,  V'  •»  H^'  =  b 1  (51) 

b  b 

2H?Y  c’  -  HTY  d'  +  hJy  b’  (52) 

0  ob  a  oa  b  ob 


Here,  w',  V',  c',  id',  and  b'  are  defined  as  are  their  counterparts 
9 

w^,  V,  c,  d,  and  b,  respectively. 

Using  (50)  and  (51),  (52)  becomes 


(Y  +  ?  )  V'  =  I’ 
a  b 


(53) 


where 
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I'  =  2H?;  Y  ,  c  .  (54) 

b  ob 

The  generalized  network  representations  of  the  apertures  correspond¬ 
ing  to  (53)  is  similar  to  that  shown  in  Figure  3,  except  for  V1  and 
I  replacing  V  and  I,  respectively.  The  reciprocity  law  then  states 


that 

-vr-*-.  ->■  t-*- 

VI'  =  V'  I  . 

This  can  be  shown  with  the  help  of  the  reciprocity  theorem. 
The  reciprocity  theorem  [2,  Section  3-8]  states  that 


(55) 


(E1  x  H2  -  E2  x  H1)  •  n  da  =  0  (56) 

J  " 

W 


where  W  is  the  closed  surface  enclosing  the  volume  containing  the 
apertures,  ( E 1 ,  H1)  and  (E2,  H2)  are  source-free  fields  within  this 
volume,  and  n  is  the  outward  unit  vector  normal  to  W.  Tet  W  be  the 
surface  consisting  of  all  the  metallic  walls  between  the  cross 
sections  of  waveguide  A  at  z  -  and  waveguide  B  at  z  -  z^  for 
some  z1  <  0  and  >  0,  and  these  two  waveguide  cross  sections,  and 
let  (E1,  H1)  and  (E2,  H2)  be  the  fields  whose  transverse  components 
are  given  by  (6)  and  (45),  respectively.  Substituting  in  (56)  and 
using  (4) ,  (56)  becomes 


i 

a 

l 

i=l 


cinaidl 


2  I 


i=l 


cXibi 


(57) 


or  in  matrix  form 

t  ->  -*T 

2d'  Y  c  =  2b  Y  ,c'  . 
oa  ob 


(58) 


The  reciprocity  law  then  follows  from  (18),  (33),  (50),  (54),  and 
(58).  The  reciprocity  law  holds  if  and  only  if  the  whole  structure 
is  reciprocal. 
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5 .  The  Scattering  Matrix 

Another  representation  of  the  apertures  is  in  terms  of  their 
scattering  matrix.  Following  Montgomery  et  al .  [10,  Section  5-14], 
the  scattering  matrix  is  defined  as 


S  = 


aa 


ba 


ab 


bb 


(59) 


where  the  iith  element  of  S  is  the  amplitude  of  the  ith  mode  in 
J  qo 

waveguide  Q  due  to  the  jth  mode  of  unit  amplitude  incident  in  wave- 
guide  0  (0  e  {A,B}).' 

Put 


a  =  [ajb  x i 

J  a 


(60) 


a'  -  [a’] 

J  b 


(61) 


It  then  follows  from  (1),  (6),  (22),  (23),  (43),  and  (45)  that 


a  «  d  -  c 


(62) 


-  c 


(63) 


The  scattering  submatrices  $aa>  ^ab*  anc^  ^bb  are  t^ien  8iven  by 
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->• 

a  =  S  c 

(64) 

aa 

-*■ 

b  =  S  c 
ba 

(65) 

(66) 

d'  =  S  c 1 
ab 

*>  ->  . 

a  =  S,,c'  . 
bb 

(67) 

The  submatrices  of  S  can  be  deduced  from  the  analysis  in 
Sections  2  and  4,  almost  immediately.  Using  (17),  (31),  and  (33), 
(62)  becomes 


a  =  (2  H  (Y  +  Y,  )  1  HTY  -  U)c  . 
a  a  b  a  oa 


(68) 


Consequently 


S  =  2  H  (Y  +  Y,)-1  HTY  -  U 
aa  a  a  b  a  oa 


(69) 


Similarly,  from  (18),  (31),  (33),  and  (65),  it  follows  that 


S,  =  2  H,(Y  +  Y,  )  1  HrY 
ba  b  a  b  a  oa 


(70) 


or,  on  using  (28)  and  (69), 


Sba-H<Saa+U)' 


(71) 


The  reciprocity  law  can  be  used  to  determine  S^.  Substi¬ 
tuting  (65)  and  (66)  into  (58),  it  becomes 


Thus 


T  T  ->  -»-T  T  -*■ 

c'L  S\Y  c  =  c  s;  Y  c* 
ab  oa  ba  ob 


-1  T 

Sah  =  Y  S7  Y  , 
at>  oa  ba  ob 


(72) 


(73) 


since  (72)  holds  for  all  c  and  c'. 
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Finally,  using  (51),  (53),  and  (54),  (63)  becomes 


['  =  (2  H,  (Y  +  1  r1  H^Y  ,  -  U)  p' 

b  a  b  b  ob 


(74) 


Consequently 


bb 


2H.(Y  +  Y,  )  1  H?  Y  . 

b  a  b  b  ob 


-  U 


(75) 


or 


Sbb 


HS  u  -  U 
ab 


(76) 


as  follows  from  (51),  (63),  (66),  and  (67). 

Let;  the  multi-mode  field  be  incident  on  the  apertures  simul¬ 
taneously  in  waveguides  A  and  B.  The  total  z-transverse  field  in  both 
waveguides  is  readily  found  from  (1)  and  (43)  a$ 


E.  = 


-Y  -2 
rai 


Y  *z 
'ai 


Y  .z 

rax 


V  c.e  d  e  ,  +  J  a,e  e  .  +  y  d!e  e  4 

f  x  — ai  £  i  —ax  L  x  -ai 

x  i  i 


■Ybiz 


if  mz 


z  <  0 


,  YYbiz 


I  bie  ^bi  +  l  cie  +  I  ai e  "  e 


•^bi 


z  >  0 


(77) 


H 


y  c.n  . 

4*  i  ai 


.z 

1  ax 


z  x  e 


ai 


-  y  a.n  . 

4  x  ax 


Y  .z 
ai 

e  z  x  e 


ai 


Y  .z 
ai 


-  I  d!n  e  ai  z  x  e  .  z  <  0 
~  i  ai  -  —ax 


-Yi  .  z  Yv  .  z 

V  i  bi  v  ?v.  bx 

>  b.a  .  e  z  x  e,  .  -  )  c  e  zxe 

j  i  bi  —  — bx  r  i  bx  — 


bi 


+ 1  e  Ybl  zx  %i  z  >  °  • 
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The  complex  power  to  the  left  of  the  apertures  is  then  given  by 


**>•  ->■  H  ->•  — ^ 

P0  Mct.td')  Yoa(c-a-d') 


(78) 


whereas  that  to  the  right  of  the  apertures  is 


-*  ->  H  *  *>■  ->  ->■ . 

Pn  —  (b  +  c  '  -fa1)  Y  u(b  -  cf  +  a  ) 
0  ob 

+ 


(79) 


Here,  as  before,  a  finite  subset  of  the  lower  order  inodes  of  cardi¬ 


nality  is  used  to  approximate  the  field  in  waveguide  Q.  By  the 
conservation  of  complex  power  law,  P^  and  P^  must  be  equal.  Put 


Y  = 


oa 


C  = 


ob 


(80) 


(81) 


Then,  using  (59),  (64),  (65),  (66),  and  (67), 


chshy*s  c  +  ch(y*s  -  s1^*)?  =  cVc  . 


(82) 


Since  C  is  completely  arbitrary,  (82)  gives 


SHY*S  +  (Y*S  -  sV)  =  Y* 


(83) 


Although  the  whole  structure  is  both  reciprocal  and  lossless. 


the  scattering  matrix  is  neither  symmetric  nor  unitary. 
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6 .  Inductive  Windows  in  a  Rectangular  Waveguide 

Consider  a  system  of  windows  located  in  the 

z=0  plane  in  a  rectangular  waveguide.  These  windows  are  assumed  uni¬ 
form  along  the  narrow  side  of  the  waveguide,  i.e.,  of  the.  inductive 
type.  The  waveguide  medium  is  assumed  characterized  by  the  real 
scalar  constitutive  parameters  c  and  y.  Figure  4  shows  the  situation 
at  hand. 

Let  a  TE-j  q  to  z  mode  of  unit  amplitude  be  incident  on  the 
windows  from  the  left.  This  mode  has  the  field  distribution 


i  -V 

E  =  e  e 


1 


i  -Yiz  a 

H*  =  ri,  e  *  z  x  e  — - —  e  rp-  x  x  e, 
—  '1  —  -1  dx  —  —1 


.i  “Y1Z 


(84) 


In  (84) 


5i  =  Ji  sin  <!  x)  i 


(85) 


and  the  subscript  "a"  is  dropped  from  and  H-,  •  Furthermore, 
it  is  assumed  that  a  <  A  <  2a  and  2b  ^ A  so  that  only  the  dominant 
mode  can  propagate  in  the  waveguide. 

Since  the  window-waveguide  structure  is  uniform  along  the 
y— axis,  and  since  the  incident  mode  has  only  an  component  that 
does  not  vary  with  y,  the  field  scattered  must  have  only  an  E^  com¬ 
ponent  that  does  not  vary  with  y.  The  only  modes  excited  in  the  wave¬ 
guide  are  therefore  TE  to  z  modes,  since  these  are  the  only  modes 

no 

having  only  an  Ey  component  that  does  not  vary  with  y  [2,  Section 
4-3] .  The  field  distribution  of  any  such  mode  is  similar  to  that 
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Figure  4.  p  inductive  windows  in  a  rectangular  waveguide. 
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of  the  dominant  mode,  viz 


E  •  e 


e 


(86) 


— n  n 


*Y_* 

H  *  +  n  e  "  zx 


x  x  e 


where 


*.  .  /Ull  V 

ab  sin  <T  *>  * 


(87) 


and 


(88) 


The  change  of  subscripts  from  i  to  n  is  for  later  convenience. 

The  total  z-transverse  field  has  only  a  y-component  of  elec¬ 
tric  field  that  does  not  vary  with  y.  It  then  follows  from  (5)  that 


each  Mm  has  only  an  x-component  that  does  not  vary  with  y: 


Mm  ■  Mm(x)x  on  l/1 


(89) 


The  generalized  network  representation  of  the  windows  can  be 
obtained  as  indicated  in  Section  2.  Here,  however,  the  whole  set  of 
modes  is  used,  which  calls  for  some  minor  changes  in  the  formulas 
th$re. 

Since  the  dominant  mode  is  the  only  incident  mode,  c  now 

becomes 


In  (90) ,  6  is  the  Kronecker  delta  function 


un 


(90) 
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6  = 


un 


u  =  n 


u  4  n 


(91) 


->■ 

The  vectors  w  ,  d,  and  b  are  likewise  vectors  of  infinite  length. 
Furthermore,  since  the  mediums  on  both  sides  of  the  windows  are 
identical , 


w  =  wv  =  w 
a  b 


(92) 


H  =  H  =  H 
a  b  o 


(93) 


Y  =  Y  ,  =  Y 
oa  ob  o 


(94) 


H  is  then  the  identity  matrix  of  infinite  order.  Consequently 


Y  =  Y  =  Y  =  H  Y  H  . 
a  b  o  o  o 


(95) 


The  generalized  network  representation  of  the  windows  then  becomes 


Y  V  =  I 


(96) 


where  Y  is  the  p  by  p  block  matrix  whose  £mth  block  is  the  matrix 


^£m  _ 


i  j J  £  m  ^  j 

J  p  xp  n=l  i 


M^e  ds1 
l  n 


W 


e  M.ds]  n 
n  i  £  m 
,m  J  p  xp 


(97) 


and  I  is  the  p  segment  vector  whose  £th  segment  is  the  vector 


1  ■  [ii]  i,  -  i-ni 

p  xl 


Mi  elds]  £ 

£  p  xl 


(98) 
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This  representation  ^s  depicted  in  Figure  5. 

Minor  changes  are  also  due  for  the  scattering  matrix  of  the 
windows.  The  higher  order  (n  >  1)  inodes  excited  are  evanescent, 
i.e.,  decay  exponentially  with  distance  from  the  windows.  Thus,  at 
sufficiently  large  distances,  only  the  dominant  mode  exists.  The 
scattering  matrix  can  then  be  defined  as 

S  = 

The  mode-amplitude  vectors  there  then  reduce  to  scalars,  viz,, 


di 


(99) 


~T  r  ->  ->  ~y  ,  .  r 

gr;  {a,  b ,  d,  a',  b  ,  d  } 


g  = 


(100) 


1  ,  g  e  (c,  c'j 


while  becomes  th^  p  segment  row  vector  h^  whose  mth  segment  is 
the  row  vector 


h™  =  [hm.]  =1-  e  Mm  ds] 

°  °Jlxpm  Jm  1  J  ’  m 


(101) 


1  X  p 


H  is  then  unity.  The  scattering  parameters  are  readily  found  from 
(69),  (71),  (73),  and  (76)  as 


m 


P  P 

-  -  -  i  -  I  [V 

1  m=l  j*l  J 


m 


e.  M™  ds 
1  J 


(102) 


b^=d,^=l+a^=l  +aj  . 


(103) 


Figure  5.  The  generalized  network  representation  of  the 


windows. 
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The  scattering  matrix  is  both  symmetric  and  unitary.  Notice 
that  the  second  term  in  the  left-hand  side  of  (82)  then  accounts  for 
the  reactive  power  basically  due  to  the  evanescent  modes. 


7 .  Capacitive  Windows  in  a  Rectangular  Waveguide 

Consider  now  the  system  of  windows  W^,  W^,  ...,  W*3  in  a  rec-T 
tangular  waveguide  shown  in  Figure  6.  These  windows,  as  can  be  seen, 
are  uniform  along  the  broad  side  of  the  waveguide,  i,e.,  of  the 
capacitive  type.  The  waveguide  medium  is  assumed  characterized  by 
the  real  scalar  constitutive  parameters  e  and  p. 

Let  a  TE^q  to  z  mode  of  unit  amplitude  be  incident  on  the 
windows  from  the  left.  This  mode  has  the  field  distribution 


i  -V 

—  *  e  *0 


i  -v 

H  =  nQ  e  z  x  eQ 


1  “V  _3_ 

jojp  £  Bx  —  x  -^0 


(104) 


where 


2  ,  ,TT  v 

e  =  l~r  sin  (—  x)  y. 
”0  /  ab  a 


(105) 


Here,  the  index  of  the  dominant  mode  is  "0"  rather  than  "1"  as  in 
(84)  and  (85).  It  is  still  assumed,  however,  that  a  <  X  <  2a  and 
2b  <  X  so  that  only  the  dominant  mode  can  propagate  in  the  waveguide. 

Since  the  window-waveguide  structure  is  uniform  along  the 
x-axis,  and  since  the  incident  mode  has  an  H  component  that  varies 
as  sin  (•-  x)  and  no  x-eomponent  of  electric  field,  so  must  be  the 
scattered  field.  The  only  modes  excited  in  the  waveguide  are 
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Figure  6.  p  capacitive  windows  in  a  rectangular  waveguide. 


36 


therefore  TE^n  to  x,  since  these  are  the  only  inodes  having  an  Hx 

TT 

component  that  varies  as  sin  (—  x)  along  the  x-axis  and  no  E  com- 

a  x 

ponent  [2,  Section  4-4].  The  z-transverse  field  of  any  such  mode 
has  the  distribution 

(106) 


±Y  z 

H  =  +  r]  e  n  zxe  +  H  y 
— n  n  —  — n  y  — 


where 


e 

— n 


e  = 
n 


n  .  /TT  v  /UTT  v 

ab~  Sln  (a  X)  COS  (T  Y)  * 


1  n  =  0 


2  n  >  1 


and 


yo 

nn  =  “  ^0 
n 


Y_  =  /  (f)2  +  (f)2  -  ^  • 


(107) 


(108) 


The  modal  expansions  in  Sections  1,  2,  and  4  are  in  terms  of 
TE  and  TM  to  z  modes.  However,  any  other  complete  set  of  modes  can 
be  utilized.  Only,  then,  the  appropriate  modal  characteristic  ad¬ 
mittances  have  to  be  used  in  (3) ,  while  everything  else  remains  un¬ 
changed.  Thus,  for  the  problem  of  capacitive  windows,  a  set  of  TE, 
to  x  modes  is  clearly  the  natural  choice,  and  (2)  and  (3)  are  to  be 
replaced  by  (108) . 
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The  total  z-transverse  field  has  only  a  y-component  of  elec- 

TT 

trie  field  that  varies  as  sin  (—  x)  along  the  x-axis.  It  then  follows 

a 

from  (5)  that  each  Mm  has  only  an  x-component  that  varies  as  sin  (~x): 

a 


Mm  =  sin  (—  x)  Mm(y)  x  on  W™  .  (109) 

—  a  — 

The  component  in  (106)  does  not  therefore  figure  in  the  analysis. 

The  generalized  network  representation  of  the  window  is  given 
by  (96)  and  depicted  in  Figure  5,  Here,  however,  Y  is  the  p  by  p 
block  matrix  whose  £mth  block  is  the  matrix 


y  r  y  -» 

1  ii  a  m 
4  P  *P 


=  t  l  \ 
n=0  i- 


sm  (—  x  )  M .  e  ds 
,  a  in 


W 


m 


e  sin  (—  x)  M^ds]  %  m 
n  a  2  P  KP 


(110) 


anfd  I  is  the  p  segment  vector  whose  2-th  segment  is  the  vector 


"  [Ii]  £  . 

p  xl 


=  [-  n. 


sin 


, IT  x 

(-x) 


M. 


e0ds) 


w 


SL  , 

p  xl 


(111) 


as  can  readily  be  verified  by  following  steps  similar  to  those  in 
Section  6. 

The  scattering  matrix  is  analogously  given  by 


S 


(112) 
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where 


r™  /7T  v  jra  , 

.  en  sin  (—  x)  M.  as 

1  0  a  i 


(113) 


m=l  j=l 


(114) 


Like  that  of  the  inductive  windows,  the  scattering  matrix  of  the 

capacitive  windows  is  both  symmetric  and  unitary.  The  second  terra 

in  the  left-hand  side  of  (82)  now  account  for  the  reactive  power  due 

to  TE.  to  x,  n  >  1,  modes. 

In 

8 .  The  Impedance  Matrix  of  the  Windows 

Let  TE^q  to  z  modes  of  arbitrary  amplitudes  c^  and  c ^  be  inci¬ 
dent  on  the  windows  (inductive  or  capacitive)  from  the  left  and  from 
the  right  respectively.  Far  from  the  windows,  only  the  same  mode  can 
exist . 


Let  v^  and  v^  be,  respectively,  the  amplitudes  of  the  E^  com¬ 


ponent  far  to  the  left  and  to  the  right  of  the  windows  referred  to 
the  z=0  plane.  It  then  follows  from  (1),  (43),  and  (99)  for  the  in¬ 
ductive  windows  and  (112)  for  the  capacitive  ones  that 


(115) 


(116) 


where  S±.  (i,j  E  {1,  2>)  is  the  ijth  element  of  the  scattering 
matrix  in  (99)  or  (112) .  The  choice  of  z=0  as  a  reference  plane  is 


only  a  matter  of  convenience.  In  matrix  form,  (115)  and  (116)  become 
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v  -  (U  +  S)  c  . 


(U7) 


In  (117) 


v  = 


V, 


(118) 


c 


Similarly,  let  i^  and  i9  be,  respectively,  the  amplitudes  of 
the  component  far  to  the  left  and  to  the  right  of  the  windows 
extrapolated  back  to  the  z=0  plane •  Then 


hoXl  (1  "  Sll)  C1  S12  c2 


^1012  S21  ci  +  C  ~  S22)  c2 


(119) 

(120) 


in  matrix  form,  (119)  and  (120)  become 


C1(J  i  «  <U  -  S)  c 


(121) 


where  is  the  characteristic  impedance  of  the  dominant  mode, 


and 


i  ss. 


"In 


(122) 


To  relate  to  network  theory,  let  (Vj  ,  -  i^)  and  ( v ^  i 2)  be 


the  complex  voltage  current  pairs  at  the  terminals  of  a  two-port 
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network  [10,  Section  4-5] .  Then 


v  =  Z  i 


(123) 


where  Z  is  the  network  impedance  matrix.  From  (117)  and  (121),  Z 
is  readily  found  as 


z  =  e1Q  (U  +  S)(U  -  S) 


-1 


(124) 


Since  S  is  symmetric,  so  is  Z.  Furthermore,  since 


z  =  c1Q(u  +  s)(u  -  s)'1  =  c10(sHs  +  S)(SHS  -  S)  1 


£10(SH  +  U)(SH  -  U)  1 


,H 


=  -  Z“  (125) 

the  elements  of  Z  are  pure  imaginary.  Finally,  using  (99),  (102), 
and  (103),  or  (112),  (113),  and  (114),  in  (124),  it  becomes 
X  X 


Z  -  3 


jX  =  -  X, 


X  X 

1  +  S11 


10  2S 


11 


(126) 


as  can  easily  be  verified  by  carrying  out  the  matrix  inversion  and 


multiplication  there. 
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9 •  The  Equivalent  Network  of  the  Windows 


The  effect  of  the  windows  on  the  dominant  waveguide  mode  is 
described  by  the  windows’  impedance  matrix  Z.  Such  a  representation 
can  be  realized  in  the  form  of  a  two-port  T-network  (10,  Section  4-5]. 
However,  because  of  (126),  the  equivalent  network  is  merely  a  shunt 
reactive  element. 

The  higher  order  modes  excited  are  evanescent,  and  are  there¬ 
fore  the  cause  for  a  localized  energy  close  to  the  windows.  In  view 
of  (3) and  (88),  the  energy  stored  close  to  the  inductive  windows  is 
predominantly  magnetic,  whereas  that  stored  close  to  the  capacitive 
windows,  in  view  of  (108),  is  predominantly  electric.  The  shunt  ele¬ 
ment  in  the  equivalent  network  is  therefore  an  inductor  for  the  induc¬ 
tive  windows,  and  a  capacitor  for  the  capacitive  windows,  as  can  be 
seen  in  Figure  7.  There 

1  + 

r.,.  - lot  inductive  windows 

10  2a1 


jX  = 


1  +  a 


310 


2ar 


for  capacitive  windows. 


(127) 


10 .  Concluding  Remarks 

In  this  chapter,  the  system  of  multiple  apertures  of  arbi- 
trary  shape  in  the  transverse  plane  between  two  cylindrical  wave¬ 
guides  has  been  considered.  The  analysis  is  based  on  the  generalized 
formulation  for  aperture  problems  [11]. 
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Figure  7.  The  equivalent  network  (a)  for  the  inductive 
windows  (b)  for  the  capacitive  windows. 
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For  a  multi-mode  incident  field,  a  representation  of  the 
system  of  apertures  in  terms  of  two  generalized  networks  in  parallel 
with  current  sources  is  obtained.  This  representation  obeys  the  two 
basic  network  laws:  the  conservation  of  complex  power  law  and  the 
reciprocity  law.  Furthermore,  each  generalized  network  depends  on 
the  modes  of  only  one  waveguide.  Thus,  for  a  given  collection  of 
apertures,  different  waveguides  can  be  considered  one  at  a  time. 

The  generalized  network  representation  can  then  be  obtained  for  any 
required  combination. 

The  scattering  matrix  of  the  apertures  is  then  deduced  from 
the  generalized  network  representation.  Although  the  aperture- 
waveguide  structure  is  both  reciprocal  and  lossless,  the  scattering 
matrix  is  neither  symmetric  nor  unitary.  This  is  because  of  the 
different  characteristic  admittances  of  the  modes,  and  the  con¬ 
sideration  of  evanescent  modes.  The  scattering  matrix  can  be  made 
symmetric  by  using  a  different  mode  normalization  from  (4),  but  no 
such  normalization  can  make  the  scattering  matrix  unitary  if  eva¬ 
nescent  modes  are  present.  A  detailed  discussion  of  this  point  is 
given  in  [12].  Since  the  scattering  submatrices  are  expressed  in 
terms  of  the  generalized  networks,  given  a  set  of  apertures,  the 
scattering  matrix  for  a  combination  of  two  waveguides  can  be  ob¬ 
tained  by  combining  the  generalized  network  for  one  waveguide  with 
that  of  the  other. 

Inductive  and  capacitive  windows  in  a  rectangular  waveguide 
are  then  considered  as  special  cases  of  the  general  problem.  The 
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windows  scattering  matrix,  since  it  involves  only  the  dominant 
waveguide  mode ,  a  propagating  mode,  is  both  symmetric  and  unitary. 
The  windows1  impedance  matrix  is  then  obtained,  apd  readily  realized 
as  a  shunt  reactive  element.  Other  problems  of  interest  that  can 
be  worked  out  as  special  cases  are  those  of  inductive  and  capaci¬ 
tive  windows  in  a  rectangular  waveguide  loaded  with  different 
dielectrics  on  both  sides  of  the  windows,  and  of  coupling  through 
small  apertures.  The  latter  is  discussed  in  a  more  general  setting 
using  related  methods  in  [13]. 

The  analysis  in  this  chapter  is  basically  theoretical,  and 
is  presented  so  that  all  the  results  and  different  relationships 
are  clearly  seen.  Results  that  relate  the  present  procedure  to 
the  mode  matching  and  conservation  of  complex  power  techniques  are 
given  in  [14]. 
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Chapter  3 


MULTIPLE  INDUCTIVE  POSTS  IN  A 


RECTANGULAR  WAVEGUIDE 


Consider  a  system  of  posts  P  ,  P 


*  •  •  • » 


PP  located  close  to 


each  other  in  a  rectangular  waveguide.  These  posts  are  assumed 
perfectly  conducting,  of  arbitrary  shape  and  thickness,  and  uniform 
along  the  narrow  side  of  the  waveguide,  i.e.,  of  the  inductive  type. 
The  medium  filling  the  waveguide  is  assumed  linear,  homogeneous, 
isotropic,  and  dissipation  free,  and  is  therefore  characterized  by 
the  real  scalar  permittivity  e  and  the  real  scalar  permeability  y. 
The  problem  considered  is  depicted  in  Figure  1. 

1.  Preliminary  Considerations 


Let  a  to  z  mode  of  unit  amplitude  be  incident  on  the 

posts  from  the  left.  This  mode  has  the  field  distribution 


E1  =  sin  (- 

y 


TT 


x)  e 


a 


x 


> 


(1) 


H 


z 


i 


where 


> 


(2) 
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Figure  1.  p  inductive  posts  in  a  rectangular  waveguide. 
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In  (2),  K  is  the  wave  number  of  the  waveguide  medium,  and  A  is  its 
wave  length.  Furthermore,  it  is  assumed  that  a  <  A  <  2a  and  2b  <  A 
so  that  only  the  dominant  mode  can  propagate  in  the  waveguide. 

Since  each  post  is  uniform  along  the  y-axis,  and  since  the 
exciting  mode  has  no  y-component  of  magnetic  field,  neither  does  the 
scattered  field.  That  is,  the  scattered  field  is  TM  to  y,  and  can 
therefore  be  derived  from  a  magnetic  vector  potential  A  having  only 
a  y-component  <f>  [2,  Section  8-7]: 


A  =  <p£  . 

The  scattered  field  is  given  in  terras  of  (p  by 


E  _  —  \J  x  V  X  (f)  y 
-  Y  2- 


H  =  V  x  (py_ 

while  <j)  itself  satisfies 

(V2  +  K2)  <j>  =  0  . 

Expanding  (4)  in  rectangular  coordinates,  the  components  of  the 
scattered  field  are  found  to  be 

4> 


(3) 


(4) 


(5) 


ES  = 
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jwe 

3y  3x 

2 
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(9  - 
v  2 
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3 

z  j(x)£  3y3z 

HS  -  3  A 

H  -  -  —  (p 
x  3z 


HS  =  0 


Hs  =  _3_  , 

h  r\  (p  , 

Z  dx 


(6) 
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Furthermore,  since  each  post  is  uniform  along  the  y-axis, 
and  since  the  exciting  mode  has  only  a  y-component  of  electric  field 
that  does  not  vary  with  y,  so  does  the  scattered  field.  It  then 
follows  from  (6)  that  (j)  is  also  independent  of  y.  The  only  components 
of  the  scattered  field  are  now  given  by 


E®  58  -  jwy  (J)  (x,z) 

Hx  "  "t  *  (XlZ) 

"z  ■  Ji  *  . 


P) 


The  total  field,  incident  plus  scattered,  must  have  zero 
tangential  electric  field  at  the  waveguide  walls.  The  incident 
field  is  a  free  waveguide  mode,  and  does  therefore  have  zero  elec¬ 
tric  field  tangent  to  the  walls.  The  scattered  field  must  then  have 
zero  tangential  electric  field  at  the  walls.  This  is  readily  accom¬ 
plished  by  setting 

4>(x,z)-0,  x  -  0,  a,  and  all  y  and  z.  (8) 

The  boundary  conditions  (8),  once  satisfied  for  any  value  of  y,  are 
clearly  satisfied  for  all  values  of  y.  Thus,  the  problem  is  basi¬ 
cally  a  two-dimensional  scalar  one  that  can  entirely  be  worked  out 
in  some  y=constant  plane  within  the  waveguide. 

In  the  next  section,  the  Green’s  function  for  the  TM  to  y 

n0 

modes  in  a  rectangular  waveguide  is  obtained.  This  is  then  used 
to  determined  <J>. 
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2. 


The  Green 1 s  Function  for  TM 

nO' 

Waveguide 


to  y  Modes  in  a  Rectangular 


Consider  a  uniform  electric  current  filament  J  directed 
across  the  waveguide  parallel  to  the  y-axis  and  located  at  (x!  ,  z1) 
as  shown  in  Figure  2. 

Since  J_  is  directed  along  the  y-axis,  the  field  produced 
must  have  only  a  y-component  of  electric  field  and  no  y-component 
of  magnetic  field.  Furthermore,  since  J  is  uniform  along  the  y- 
axis,  so  must  be  E^.  This  can  readily  be  established,  for  instance, 
by  the  reciprocity  theorem  [2,  Section  3-8],  A  magnetic  vector 
potential  having  only  a  y-component  <j),  proportional  to  E^,  as  is 
seen  in  Section  1,  can  then  be  used  to  derive  all  field  components. 

Only  TMnQ  to  y  (TE^q  to  z)  modes  can  be  excited  in  the 
waveguide,  since  these  are  the  only  modes  having  only  an  E^  com¬ 
ponent  that  does  not  vary  with  y  and  no  component  [2,  Sections 
3-4, 4].  The  potential  function  ({)  due  to  the  filament,  relabeled  G, 
is  therefore  referred  to  as  the  Green1 s  function  for  TMnQ  to  y 
(TE^q  to  z)  modes  in  a  rectangular  waveguide.  Below,  G  is  found  as 
a  series  of  these  modes. 

The  wave  equation  satisfied  by  G,  for  each  y,  is 


(-^~2  +  +  K2)  G(x,z)  -  -  6 (x-x 1 )  6(z-zf)  (9) 

3x  3z 

subject  to  the  boundary  conditions 


=  0,a , 


G(x,z)  =  0, 


x 


and  all  z. 


(10) 


50 


Figure  2 


tx 


— ►  2 


An  electric  current  filament  J  in  a  recgangular 
waveguide  parallel  to  the  y-axis. 
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In  (10),  6  is  the  Dirac  delta  function.  Multiplying  throughout 

(9)  by  sin  (—  x) ,  then  integrating  over  x  from  0  to  a,  it  becomes 
a. 


(M)2) 

a 


a 

r 


0 


G(x,z) 


sin 


dx 


=  -  sin  (: —  x!)  6(z-z!)  .  (IX) 

a 


Put 


Gn(Z) 


a 

G(x,z)  sin  (—  x)  dx 
3 

0 


n  =  1,2, 


(12) 


The  one-dimensional  wave  equation  (11)  then  becomes 


Y  )  G  (z) 
n  n 


sin  ( —  x')  6(z-z')  . 
a. 


(13) 


For  the  solution  of  (13)  to  represent  waves  traveling  away 


from  the  filament,  G  must  be  of  the  form 

n 

f  -Y  Z 

A  n  .  . 

A  e  z  >  z 

n 


G  (z) 
n 


< 


B 

n 


z  <  z 1 


(14) 


where  A  and  B  are  constants  to  be  determined.  Since  G  is  pro- 
n  n  * 

portional  to  E  ,  it  is  continuous  across  the  filament  at  z=z! 

[2,  Section  1-14],  and  so  is  G  .  Thus 

n 


-Y  z ' 
n 


Y  Z  1 

*n 


A  e 
n 


B  e 
n 


=  0  . 


(15) 
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Furthermore,  integrating  (13)  over  z  from  z'-A  to  z'+A,  then  let¬ 
ting  A  go  to  zero,  it  becomes 


£  G„<‘> 


+ 
,  I 


,m  f  N 

SID  ( -  XT) 

a 


(16) 


That  is,  —  G  is  discontinuous  at  z  =  z*  by  the  amount  -  sin  (—  x1). 
dz  n  a 


Thus 


Ay  e  ^  +  B  y  en*  =  sin  (—  x')  . 


n  n 


n  n 


(17) 


Solving  (15)  and  (17)  simultaneously,  A  and  B  are  found 

n  n 


to  be 


An=  2TS±n  (f  X’}  e'n 
n 


Y  z ' 


(18) 


-y  z ' 


r,  1  •  /Mi  n  '  n 

Bn  W  S1"  (T  *  >  6 

n 


(19) 


Combining  (18)  and  (19)  with  (14),  G  becones 

n 


-y  | z-z 1 1 

G  (z)  =  — —  sin  (—  x1 )  e  n  ,  n«l,2f ,  .  . 

n  zy  a 

n 


(20) 


By  Fourier  theory  [15,  Section  43  ]  ,  (12)  can  be  inverted  as 

-yJz-z'  | 


g(x,z i x' ,z ' )  =  —  y 


.  /nTT  v  .  /nTT  'n* 

s  in  ( —  x  )  s  m  ( —  x  )  e 
a  a 


n=l 


(21) 


Clearly,  G  satisfies  the  boundary  conditions  (10). 
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3,  Basic  Formulation 

Let  (E*,  H^)  be  incident  while  all  the  posts  are  absent, 
and  (_E (J)  >  H(J) )  be  the  field  produced  by  an  electric  current 

of  density  _J  «  Jm,  where  J™  is  the  current  on  Pm,  while  all  the 
m=l 

posts  are  absent.  By  the  uniqueness  theorem  [2,  Section  3-3], 

(E*  + _E(J) ,  H*  +  H(J))  is  identical  with  the  original  field 
whenever 

nm  x  (E1  +  E(J))  =0  on  P®  .  (22) 

In  (.22),  is  the  outward  unit  vector  normal  to  P™.  (E(J)»  HQ)) 

must  then  have  the  field  distribution  (7).  Since 


m 


n  x  (H  (J ) 


H  (  J) 


v=o. 


)  =  J  on  P 


m 


(23) 


v=0 


where  V  is  the  distance  along  n.m  from  Pm,  d™  has  only  a  y-component 
that  does  not  vary  with  y: 


Tm  \ 

J  -  J  (x,z)  y_  . 


(24) 


As  is  pointed  out  in  Section  1,  the  problem  is  a  two- 
dimensional  scalar  one  that  can  be  worked  out  in  some  y=constant 
plane  within  the  waveguide.  Thus,  all  source  and  field  points  are, 
hereafter,  assumed  located  in  any  such  plane. 

By  definition,  G (x, z | x r , z 1 )%  is  the  magnetic  vector  poten¬ 
tial  produced  at  any  point  (x,z)  by  a  unit  electric  current  fila¬ 
ment  in  the  y-direction  located  at  (x’^z*).  By  superposition,  then, 
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v[  produces  at  (x,z)  the  magnetic  vector  potential  <J>  where 


4>(x,z)  = 


Jm(x',z')  G(x,z|x',z')  dt 


m=l  ^,m 


(25) 


dt' 


/ 


(dx')2  +  (dz’)2 


In  (25),  G  is  given  by  (21),  and  primed  and  unprimed  coordinates 
denote,  respectively,  source  and  field  points. 

Since  G(x,z|x,,zT)  is  a  solution  of  the  homogeneous  wave 
equation  (5)  for  all  (x,z)  ^  (x^z1),  so  is  (j).  Furthermore,  cj)  satis¬ 
fies  the  boundary  conditions  (8)  by  virtue  of  (10).  Thus  (f>,  and 
consequently  the  complete  field  solution,  can  be  found  once  all  Jm 
are  known.  Using  (1),  (7),  and  (25),  (22)  becomes 


2  p  r 

-  jcoy  Jm(x’,zl)  G(x,z|xr,zf  dtT  «  0, 

m=l  Jm 


,TT  \ 

sin  (—  x)  e 
a 


Q 

(x,z)  e  c  ,  1  £  i  <_  p 


(26) 


which  is  an  integral  equation  for  J_. 

The  higher  order  (n  >  1)  modes  are  evanescent,  i.e.,  decay 
exponentially  with  distance  from  the  posts.  Thus,  at  sufficiently 
large  distances,  only  the  dominant  (n  =  1)  mode  can  exist  in  the 
waveguide.  The  reflection  coefficient  of  the  dominant  mode  is 
readily  found  from  (7),  (21),  and  (25)  as 


•  Z77  I  \ 

sm  (—  x  )  e 
a 


r 


1  m=l  Jm 

Vj 


dt'  . 


(27) 
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The  transmission  coefficient  of  the  dominant  mode  is  then 


T  =  1 


_ 


ay 


1  m=l 


Jm(x 1 ,  z  1 )  sin 


m 


/7F  i  \ 

(—  xf ) 
a 


Yj_z 


dt* 


(28) 


4.  The  Scattering  Matrix 

Following  Montgomery  et  al.  [10,  Section  5-14],  the  scatter¬ 
ing  matrix  of  the  posts  is  defined  as 


S  = 


11 


21 


12 


22 


(29) 


In  (29),  S^  and  S^  are,  respectively,  the  amplitudes  of  the  dominant 
mode  reflected  to  the  left  and  transmitted  to  the  right  of  the  posts 
due  to  an  incident  TE^  to  2  mode  of  unit  amplitude  from  the  left.  Con¬ 
sequently,  S^  and  are  given  by  (27)  and  (28),  respectively. 

Similarly,  S^  and  S^  are,  respectively,  the  reflection 
and  transmission  coefficients  of  a  TE^q  to  z  mode  of  unit  amplitude 
incident  on  the  posts  from  the  right.  This  mode  has  the  field  dis¬ 
tribution 


E  = 

y 

E1  = 

X 


H  = 


Yi  2 

/TT  x  f  1 
sm  (—  x)  e 
a 


Y1  TT 

- —  sin  (—  x)  e 
jwy  a 


nz 


-TT  ,71  x 

- -  cos  (—  x)  e 

joopa  a 


nz 


(30) 
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The  previous  analysis  carries  through  in  this  case.  Thus, 
the  scattered  field  is  given  by  (7)  and  (25),  but  with 

_JT  =  J^rm,  now  replacing  _J  in  (25),  determined  by  solving  the 


m-1 

integral  equation 


TT  Y1Z 

sin  (—  x)  e  -  jo)jj 
a 


m=l 


J|IT1(x,,zt)  G(x,z|x!,zr)  dt 1  =  0, 


m 


£ 

(x,z)  £  C  ,  1  <  £  <  p 


(31) 


rather  than  (26).  It  then  follows  from  (7),  (21),  and  (25)  that 


_  2m 


22  ay 


1  m=l  l 


m 


v  z, 1 

-j.  •  m »  .  .  v  .  ,77  *  v  1 1  ,  , 

J  (x  ,z  )  sm  (—  x  )  e  dt 

a 


(32) 


ay 


1  m=l 


m 


-y  z! 

J,m(xT,zT)  sin  (“  x 1 )  e  dtf  .  (33) 

a 


is 


The  scattering  matrix  is  both  symmetric  and  unitary.  That 


S  = 


ssH  =  sHs  =  u 


(34) 


where  T  and  H  denote,  respectively,  matrix  transpose  and  Hermitian, 

and  U  is  the  identity  matrix. 

11  2  2 

Let  (E  ,  H.  )  and  (E  ,  H  )  be  the  z~transverse  fields  in  the 
waveguide,  sufficiently  far  from  the  posts,  due  to  to  z  modes 

of  arbitrary  amplitudes  c^  and  c^  incident  from  the  left  and  from 
the  right  of  the  posts,  respectively.  It  then  follows  from  (1),  (7), 
(21),  (25),  (27),  (28),  (30),  (32),  and  (33)  that 
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-v 


V 


E  =  < 


c  (e  +  Sn ..  e  )  sin  (—  x)  y  z  «  0 

l  il  a 


-Y-,  z 

TT  1 

C1S21  Sin  6  —  z  »  0 


H  = 


E  = 


"V  Y1Z  77 

-TL  c.  (e  -  S.,  e  )  sin  (—  x)  x  z  «  0 
11  11  a  — 


-Y.  Z 

„  .TT  .  '1 

-n,  c  S  sm  (—  x)  e  x 
1  1  21  a  — 


YjZ 

C2S12  Sin  (5  X)  e  * 


z  »  0 


z  «  0 


(35) 


V 


-YiZ 


c2 (e  +  S22  e  )  sin  (—  x)  y  z  »  0 


Y  -t  Z 

„  .IT  v  '1 

n,  c  S  sm  (-  x)  e  x 

121/  a  — 


z  «  0 


H  = 


V 


-v 


n1c  (e  -  S22  e  )  sin  (—  x)  x  z  »  0. 


In  (35)  ,  is  the  characteristic  admittance  of  the  dominant  waveguide 
mode: 


(36) 


Let  W  be  the  closed  surface  consisting  of  all  metallic  walls 
between  the  two  waveguide  cross  sections  at  z  =  z^  and  z^>  for  some 
z^  <<  0  and  »  0>  and  these  two  cross  sections.  The  reciprocity 
theorem  then  states  that 


r  i  ?2i 

( E  x  H-  !  X  H  )  •  n  ds  =  0 

w 


(37) 
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where  n  is  the  outward  unit  vector  normal  to  W.  Substituting  (35) 
into  (37),  there  then  results 


abn1c1S12c2  abn1c1S21c2 


(38) 


whence 


S12  "  S21  ’ 


(39) 


The  scattering  matrix  is  symmetric  if  and  only  if  the  whole  structure 
is  reciprocal. 

That  S  is  unitary  follows  from  conservation  of  power  con¬ 
siderations.  Let  the  two  dominant  modes  be  simultaneously  incident 
on  the  posts  from  the  left  and  from  the  right.  The  complex  power 
scattered  far  to  the  left  and  to  the  right  of  the  posts  is  basically 


P 

sc 


ab 


n1(l 


C1S11  +  C2S1 2 ' 


+ 


'  C1S21  +  C2S22 


I2) 


(40) 


whereas  that  incident  is  given  by 


(41) 


Since  the  structure  is  lossless,  and  since  Pin  and  Psc  are  real, 
they  must  be  equal .  Put 

C1 


(42) 


Then 


ab  -*H+ 

T  ni  c  c 


+H 

c 


SHS  c 


(43) 


or 


sHs  =  U. 


(44) 
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5.  The  Impedance  Matrix 

Let  TE^q  to  z  modes  of  arbitrary  amplitudes  and  c 2  be 
incident  on  the  posts  from  the  left  and  from  the  right,  respectively. 

Let  and  v 2  be,  respectively,  the  amplitudes  of  the  com¬ 
ponent  far  to  the  left  and  to  the  right  of  the  posts  referred  to  the 
z  =  0  plane.  It  then  follows  from  (35)  that 


(1  +  S1;L)  c±  +  S12  c2 

(45) 

=  S21  cx  +  (1  +  S22)  c2  . 

(46) 

The  choice  of  z  =  0  as  a  reference  plane  is  only  a  matter  of  con¬ 
venience.  In  matrix  form,  (45)  and  (46)  become 


v  =  (U  +  S)  c 


(47) 


where 

(48) 

Similarly,  let  i-^  and  i 2  be,  respectively,  the  amplitudes 
of  the  component  far  to  the  left  and  to  the  right  of  the  posts 
extrapolated  back  to  the  z  =  0  plane.  Then 


'^l1!  (1  Sll)  C1  ~  S12C2 


^2  =  “  +  (1  “  ° 


21  1 


22  2 


(49) 

(50) 


In  matrix  form,  (49)  and  (50)  become 


i  -  (U  -  S)  c 


(51) 
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where 


(52) 


To  relate  to  network  theory,  let  (v^,  -  i^)  and  (v^,  i^) 
be  the  complex  voltage-current  pairs  at  the  terminals  of  a  two-part 
network  [10,  Section  5-2].  Then 


v  =  Z  i  (53) 

where  Z  is  the  network  impedance  matrix.  From  (47)  and  (51),  Z 
is  readily  found  as 


Z  =  q(U  +  S)(U  -  S) 


(54) 


Since  S  is  symmetric,  so  is  Z.  Furthermore,  since 
Z  =  cx(u  +  S)(U  -  s)'1  =  q  (SHS  +  S)(SHS  -  S)_1 

-  (SH  +  U)(SH  -  U)"1 

=  -  zH 

the  elements  of  Z  are  pure  imaginary.  Thus 
X11  X12 

Z  =  j 

X21  X22 

—  J 

X12  "  X21 ’ 


(55) 


(56) 
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6 .  The  Equivalent  Network 

The  complete  field  solution  is  seldom  needed.  Rather,  the 
effect  of  the  posts  on  the  dominant  waveguide  mode  is  what  must 
accurately  be  described.  From  an  engineering  perspective,  a  de¬ 
scription  in  terms  of  a  network  of  lumped  elements  is  preferred. 

The  effect  of  the  posts  on  the  dominant  waveguide  mode  is 
fully  described  by  the  posts’  impedance  matrix  Z.  Such  a  representa¬ 
tion  can  be  realized  in  the  form  of  a  two-port  T-network  [10,  Sec¬ 
tion  4-5]. 

The  characteristic  impedances  of  the  TM^  t0  y  modes  are 

given  by 


(57) 

Since  these  modes  are  evanescent,  the  energy  stored  close  to  the 
posts,  in  view  of  (57),  is  predominantly  magnetic.  This  effect  can 
suitably  be  represented  by  an  inductor  in  the  shunt  arm  of  the  net¬ 
work.  The  elements  in  the  series  arms,  however,  are  capacitors  to 
account  for  the  charge  difference  across  the  posts  in  the  z  direction. 
The  equivalent  network  of  the  posts  is  shown  in  Figure  3. 
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H^if  X|^)  j  t  ^22~  ^12^ 


Figure  3.  The  equivalent  network  of  the  posts. 
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7 .  Solution  of  the  Integral  Equation 


The  integral  equation  (26)  can  be  put  in  the  compact  form 


Zm(Jra)  -  V 


m=l 


_m/Tnk  , 

Z  (J  )  =  jwy 


Jm(x' 9z')  G(x>z\x' 9z')  dtT 


m 


V  =  sin  (—  x)  e 
a 


’V 


£ 

(x,z)eC,l<£<p. 


(58) 


An  exact  solution  of  (58)  can  rarely  be  obtained,  and  an  approximate 

solution  has  then  to  be  sought. 

Let  each  C™  be  approximated  by  a  polygon  £  m  =  {  S™,  S™  >  •  •  •  j  S^} 

qm 

as  shown  in  Figure  4,  and  put 


m 


Jm(x' ,z') 


In  (59),  I?  are  complex  coefficients  to  be  determined,  whereas  each 

m  ^  m 

J.  is  a  real  function  that  vanishes  on  all  S.  ,  .  ,  but  is  otherwise 
3  1  t  3 

unspecified.  Substituting  (59)  into  (58),  it  becomes 


\  l“  j?(x’,z')  . 
j=l  2  2 


(59) 


m 

?  \  im  zm(Jm 

mti  j=l  1  J  J 


)  +  r  =  V  , 


(x,z)  e  sL  1  <  i  <  p,  1  <  i  <  q£ 


(60) 


where  Z™(J™)  is  given  by  (48)  except  that  the  integration  is  taken 
over  S™,  and  r  is  a  residual  term.  A  Galerkin  solution  [4,  Sec¬ 
tion  1-3]  can  be  obtained  by  requiring  that  r  be  orthogonal  to 

all  Jm  . 

J 
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Figure  4. 


Cm  approximated  by  a  polygon  \ 


m 
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Define  the  inner  product 


<A,  B>  = 


A  B  dt  . 


(61) 


(!)  t 

i-1 


Taking  the  inner  product  of  (60)  with  each  »  and  enforcing  the 
Galerkin  condition 


<r ,  J±>  =  0  , 


1  <  £  <  p  ,  l<i<q 


(62) 


there  then  results  the  system  of  equations 


m 


,.m  ,„m,  .in.  _Jl.  T 
I .  <Z. (JJ  ,  j1>  =  <V,  Ji>  , 


m=l  j=l  J  J  -1 


1  <  £  <  p  ,  l<i<q 


(63) 


In  matrix  form,  (63)  becomes 


=  ->-  ->■ 
Z  I  =  V 


(64) 


Vhere  Z  is  a  p  by  p  block  matrix  whose  £mth  block  is  the  matrix 

2  '  1Zij]  l  m‘  I<z3<Jj)-  V1  1  m  (65) 

J  q  xq  j  j  q  xq 

and  I  and  V  are  p  segment  vectors  whose  mth  and  £th  segments  are 
the  vectors 


i"  -  [i"i 

]  ,"xl 


(66) 


V2  -  [V2]  ,  -  [<V,  J2>] 

q  xl  q  xl 


(67) 


respectively. 
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The  currents  Jm  given  by  (59),  with  the  coefficients  I? 

3 

determined  from  (65),  form  the  Galerkin  solution  of  (58).  A  Galerkin 
solution  of  (31)  can  be  obtained  in  a  similar  manner.  Clearly,  then, 
using  the  same  J™,  the  solution  is  given  by  (59),  but  with  the  coef¬ 
ficients  now  determined  by  solving  (64)  with  the  right-hand  side 

— y 

vector  V  conjugated. 


8 .  Evaluation  of  the  System  of  Equations 

The  construction  of  Z  in  (64)  constitutes  a  large  portion 
of  the  work  involved  in  the  numerical  solution.  An  efficient  evalu¬ 
ation  of  the  elements  of  Z  is  therefore  necessary  for  the  success  of 
the  solution. 

The  ij th  element  of  the  £mth  block  of  Z  is  given  by 


Jim 

z±.  =  jo>y 


f 


J^(x,z)  dt  I  J^x’jZ1)  G(x,z|x,,zT)  dt 1  (68) 


sm 


where  J™  are  so  far  unspecified.  A  particularly  simple  choice  for 

J .  is 
3 


Jj(x\z’)  = 


1  (X’.z')  e  S" 


(69) 


0  (x 1 ,  z  1 )  e  S  . 


i  *  j 


m 

which  corresponds  to  a  pulse  expansion  of  J  .  ^±3  ^^ien  becomes 


_&m 

z. .  =  3 wy 

13 


f  dtf 

s£  s? 

1  3 


G(x,  z  |x ' , z '  )  dt '  . 


(70) 
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Put 


£m,  . 

(x,z) 


G(x ,  z  |x ' ,  z  '  )  dt ' , 


sm 

J 


(x,z)  e 


(71) 


By  the  first  mean  value  theorem  of  integration  [16,  Section  7-18], 

i 

there  exists  a  point  (x  ,z  )  e  S.  such  that 

o  o  1 


„£m  .  _  £  „£m,  . 

Z..  =  jcoy  L.  R  .  (x  ,z  ) 
1J  J  1  ij  o’  o' 


(72) 


where 


T  £  .  £  £.2  ,  £  £.2 

1i "  (xi+i  ‘  V  +  (Vi  -  2i> 


(73) 


is  the  length  of  S^. 


,£m 


The  evaluation  of  Z  is  now  completed  by  integrating  G 


„m 


over  S . .  Put 
1 


G(x  , z  x',z')  =  (G  +  G  +  G  )(x  , z  x',z') 
o  o  0  1  zoo 


(74) 


In  (74) 


r  (  I  i  in  1  .  ,TT  .  .  (T[  ,  a  IVZ  ^1 

G0(Vzo|x  *Z  }  =T^Sin  (7Xo)  Sin  (a  x  )e 

00  _  HU  _z. 

G1  (x  ,z  |x'  ,z’)  =  ~~  £  —  sin  (—  x  )  sin  (—  x')e  a  ° 

1  v  o  o  1  7T  . ..  n  a  o  a 

n=l 

_  JIl  Z  -z  T  I 

G0(x  ,  z  |xf  ,z?)  =  —  [^sin  (—  x  )  sin(~  x!)  e  a  ° 

2  o  o 1  tt  a  o  a 


r  .  /Till  x  .  /IITT  ,  N 

4-  )  sm  ( —  x  )  sm  ( —  xT) 

L  a  o  a 

n=2 


x(- 


-  ^|z  -z'lB 

a  o  n 


a.,z  .z. 

a  I  o 


n 


-)i 


(75) 
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.  /  2  ,ir.  2  .  tt  /,2a.  2  .  TT 

Yl  -  J  /  <  -  (-)  '  3  -  /  (T)  -  1  *  3  I  h 


,nrrN2  2  it  /  2  ,2a.  2  tt  Q  _ 

(— )  ~  K  =  ~  /  n  -  (-r-)  >  n  >  2 

a  a  /  A  an  — 


(76) 


The  decomposition  (74),  therefore,  amounts  to  expressing  the  dynamic 
Green’s  function  G  in  terms  of  a  dominant  mode  wave  G^,  the  cor¬ 
responding  static  Green’s  function  G^ ,  which  can  be  obtained  from 
(21)  by  setting  K  equal  to  zero,  plus  correction  terms  G^. 

The  series  defining  G^  is  readily  summed  to  give 

.j  cosh (“|  z  — z  ’  |  )  -  cos  —  (x  +x’) 

G1  (x  , z  |x 1 , z ' )  *7- r  log  ( - = — ° - 1 - - - )  (77) 

cosh(—  z  -z ' | )  -  cos—  (x  -x') 
a  o  a  o 


where  log  denotes  the  natural  logarithm.  The  details  of  the  summa¬ 
tion  are  given  in  Appendix  A.  The  series  in  G^  is  dominated  by  an 
exponentially  convergent  series  of  positive  monotonically  decreasing 
terms  (see  Appendix  B) ,  and  can  therefore  be  summed  directly  at  a 
minimal  cost. 

The  integration  of  G  can  be  carried  out  numerically,  and 
for  that  purpose  any  quadrature  rule  can  be  used.  Thus 


; 


sm 

1 


G(x  ,z  |xf ,z’)  dt ’ 
o  o 1 


Lm  N 

-fl  ’k<G0  +  Gl+G2)(Vzo! 

k=l 


m  /1  \m.  m 

(1  -  pt)xj  +  PkVl.  (1  -  Pk)Zj  +  PkZ.+1 


(78) 
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In  (78),  N  is  tf*e  order  of  the  rule,  are  its  coefficients, 


and  p,  determine  the  location  of  its  abscissas. 


When  evaluating  the  diagonal  elements  of  Z,  Z^,  ^  °^^ers 

l 

a  logarithmic  singularity  at  (x^jZ^)  e  S_^  that  requires  particular 
attention.  In  Appendix  A,  the  singular  part  of  is  found  to  be 


Gls(V2olX’’Z,) 


1  log  (J  / (x'~x0)2  t  (z'-zQ)2) .  (79) 


2tt 


Put 


G  (x  ,z  |x',z’) 
lp  o  o' 


(80) 


Then 


G(x  ,z  |x',z')  dt' 
o  o 1 


sf 

l 


G  (x  ,z  |x',z')  dt'  + 
Is  o  o' 


(Gq"K^1P+G2)  (xo » z0  lx ' 1  z '  )dt ' 


■  -  s  [L  l0*(rr-7> +  Li(los  <L1  -  L»  - 

L  #  —  Li 
1 


+ 


dt’  . 


(81) 


i 

£  £ 

Here,  L  is  the  distance  between  (x^9z^)  and  (x^,  z^) .  The  integral 
on  the  right-hand  side  of  (81)  has  no  singularity  at  (x0>z0)*  and 
can  therefore  be  evaluated  through  (78). 
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The  ith  element  of  the  £th  segment  of  V  is  given  by 


£  f  £  ^1Z 

V.  =  J.(x,z)  sin  (—  x)e  dt 

1  Z  1  a 

s. 

x 


(82) 


which,  upon  using  (69),  becomes 


1 


s£ 

1 


rTT  ,  ^1  dt. 

sin  (—  x)  e 

a 


(83) 


The  integration  in  (83)  can  be  carried  out  exactly.  However,  a 

point  (x  ,z  )  exists  such  that 
o  o 


V .  =  L ,  sin  (—  x  )  e 
1  x  a  o 


-Vo 


(84) 


Actually,  finding  such  points  (xq,zo)  and  (xq,zo)  is  at  least 
as  difficult  as  computing  the  integrals  themselves.  For  sufficiently 
small  L^,  however,  the  mid-point  of  S ±  can  replace  these  points  while 
introducing  only  very  little  error.  The  system  of  equations  thus 
obtained,  clearly,  is  one  that  results  from  enforcing  the  point  match¬ 
ing  condition 


r(x,z)  =  0, 


£  ,  £  £  ,  £ 

X  .  .  ,  +  X  .  z  1  +  z  0 

(x,z)  e  {(-^2 - “  2 - )l  1  1  1  1  P.  1  1  i  1  qh 

(85) 


in  (60)  rather  than  the  Galerkin  condition  (62) 


9 .  Numerical  Results 

The  solution  procedure  presented  is  readily  translated  into  a 
computer  program.  The  elements  of  the  scattering  matrix  and  the 
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reactances  of  the  equivalent  T-network  are  basically  the  parameters 
to  be  computed. 

The  scattering  parameters,  thanks  to  (59)  and  (69),  are  com¬ 
puted  by 


11 


21 


22 


12 


jwH 

ay. 


m 

I  \ 

1  m=l  j=l 


t?  i 
J  J 


sm 

3 


-y  1 

sin  (—  xf)  e  ^  dt? 
a 


m 

=  1  -  Jai  1  \  I 

aY, 


1  m=l  j=l 


sm 

J 


TT  Y1Z' 

sin  (—  xT)  e  dt1 
a 


m 


_  Jml  x  \  I«m 

aYl  m=l  J-l  j  l 


m 

J3 


Yt  z 

.IT  ,s  1  , 

sin  (—  x  )  e  dt 
a 


m 


.m 


!  _  iMi  y  1  i: 
aYl  m=l  A  J  J 


sm 

J 


-y-  z’ 

sin  (—  x 1 )  e  dt  * 

a 


(86) 


where  I1!1  and  l!m  are  the  solutions  of  (64)  with  right-hand  side 

J  1 

vectors  V  and  V  ,  respectively.  The  impedance  matrix  is  then  com¬ 
puted  through  (54) .  In  carrying  out  the  integrations  in  (86) ,  and 
also  in  (78)  and  (81),  an  eight-point  Gauss-Radau  quadrature  rule 
[17]  is  used.  Table  1  shows  the  p^  and  q^  of  the  rule. 

Because  of  the  approximations  involved  in  the  solution,  however, 
the  scattering  matrix  need  no  longer  be  symmetric  nor  unitary.  To 
determine  the  impedance  matrix,  S.^  and  S21  are  first  replaced  by 
their  average 


S 

av 


+  s21) 


(87) 


The  impedance  matrix  can  then  have  a  non-zero  real  part. 
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To  test  the  solution,  the  computer  program  is  run  for  a  few 
selected  problems.  In  particular,  the  problems  of  the  circular  post, 
of  the  symmetrical  thin  window,  and  of  the  triple  circular  post  are 
considered.  Some  of  the  results  obtained  are  plotted  in  Figures  5-9. 

In  all  the  cases,  the  convergence  for  the  inductive  reactance 
is  monotonic  and  from  above,  as  can  readily  be  seen  from  Figures  6 
and  8.  The  computed  reactances  are  found  to  agree  well  with  the  data 
in  the  Waveguide  Handbook  (WGHB )  [3],  with  only  a  few  segments  needed 
even  for  large  posts.  A  complete  assessment  of  the  solution  perform¬ 
ance  should  also  consider  the  (Frobenius)  norm  of  the  real  part  of  Z 
and  the  modulus  of  difference  in  transmission  coefficients.  These 

—8 

two  numbers  are  computed  in  all  program  runs,  and  are  usually  0(10  ). 

Perhaps  the  most  interesting  observation  can  be  drawn  by 

examining  Figure  5  for  the  centered  circular  post.  For  large  posts 
d  \X1 

(—  >  0.25),  —  -r—  is  no  longer  frequency  independent  as  is  the  case 

a  £-j_  /a 

with  smaller  posts  (^  <_  0.25),  but  rather  branches  out.  Figure  9 
for  the  symmetrical  triple  circular  post  displays  yet  another  almost 
frequency  independent  characteristic.  This  is  not  surprising,  however, 
since  this  configuration  cancels  out  the  first  six  higher  order  modes 
[18,  Section  5-1.3]  . 

10.  Concluding  Remarks 

The  system  of  inductive  posts  in  a  rectangular  waveguide, 
i.e.,  of  all  the  metallic  obstacles  that  are  uniform  along  the  narrow 
side  of  the  waveguide,  but  are  otherwise  of  arbitrary  shape 
and  thickness,  has  been  considered  in  this  chapter. 
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Figure  5.  Network  reactances  of  the  centered  circular  post.  The 
number  of  segments  used  is  24. 
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Figure  6.  The  convergence  behavior  for  the  centered  circular 

post  (—  =  1.2,  —  =  0.1). 
a  a 
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x,  Xi 


Figure  7.  Network  reactance  of  the  symmetrical  thin  window.  The 
number  of  segments  per  diaphragm  is  24. 


py  1° 
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Figure  8.  The  convergence  behavior  for  the  symmetrical  thin 

windows  (y  =  0.8,  —  =  0.5). 

A  a 


n 


Figure  9.  Network  reactances  of  the  symmetrical  triple  post.  The 
number  of  segments  used  for  each  post  is  16, 
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A  complete  field  analysis  of  the  problem  i s  first  given.  From 
this  analysis,  the  scattering  and  impedance  matrix  representations  of 
the  system  of  posts  fully  describing  its  effect  on  the  dominant  wave^ 
guide  mode  are  obtained.  Since  the  whole  structure  is  both  recipro¬ 
cal  and  lossless,  the  scattering  and  impedance  matrices  are  symmetric 
and  unitary,  and  symmetric  and  pure  imaginary,  respectively.  The 
latter  is  then  real  ized  in  the  form  of  a  T-network  of  reactive  ele¬ 
ments  , 

The  reactances  of  some  post  configurations  are  computed  in 
Section  9,  and  more  results  can  be  found  in  fl9|.  In  the  actual 
computation,  pulse  expansions  of  the  currents  induced  on  the  posts 
are  used.  Although  chosen  primarily  so  as  to  render  the  procedure 
economical,  the  choice  is  very  natural,  since  pulses  are  instru¬ 
mental  in  the  definition  of  integration  [16,  Chapter  10].  This 
choice  has  proven  very  successful,  nevertheless,  as  is  evident  by 
the  performance  of  the  solution. 

The  circular  post  problem  was  first  treated  using  the 
Variational  Method  f 1 ,  Chapter  2],  [2,  Section  8-7].  The  Vari¬ 
ational  Method  was  also  applied  to  solve  the  problems  of  the 
inductive  thin  and  thick  irises  |9,  Sections  8-4  and  8-5],  and 
of  an  array  of  symmetric  thick  irises  [23|,  Single  and  triple 
circular  posts  have  been  considered  in  (7],  I  8],  where  many  of 
their  characteristics  have  been  discussed.  The  present  analysis, 
however,  is  quite  general.  It  can  also  become  the  first  step 
in  the  solution  of  a  system  of  dielectric  posts  in  the  inductive 
position  in  a  rectangular  waveguide. 
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Appendix  A 

Consider  the  function  defined  by  the  series 


G1  "  I 

n=l 


sin(nri)  sinCnti'je  n° 


a  =  k  -  V 


(A.l) 


Since 


sin(nn)  sin(nri')  =  77  (cos  n(n~n’)  ”  cos  n(n+n')) 


A  Re(e~^n(n-r| ' }  -  e-3n<rK-n')) 


(A. 2) 


(A.l)  becomes 


G1  ■  2  Re(  l,  n 
n=l 


1  e-n(a+j (n-n ' ) ) 


uu 

y  1  e-n(o+j (n+n')) 


).  (A. 3) 


n=l 


In  (A. 2)  and  (A. 3),  Re(z)  denotes  the  real  part  of  z. 
Since 


oo 

=  I  zn»  Z  e  C,  |z|  <  1 


n=0 


(A. 4) 


and  the  series  in  (A. 4)  converges  uniformly  for  all  z,  |  z  |  _<  |  z  |  <  1, 
a  term  by  term  integration  can  be  carried  out  [20,  Section  3^4], 
giving 


z 

f° 


dz  _  y 

1-z  L 
n=0 


n 

z  dz 


n+1  n 

oo  2  oo  ^ 

-  log  (1-z  )  -  I  -2-  -  1-2. 

°  n=0  n+1  n-1  n 


(A. 5) 
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In  (A. 5),  log  denotes  the  natural 

-(a+j  (n-n ’ ) )  .  -(cf+j  (n+n')) 

e  and  e 

in  (A. 3),  becomes 


logarith 
in  (A. 5), 


Putting  z^  equal  to 
then  using  the  results 


G1  =  y  Re [log 


j  _  e~(o+j (n+n')) 
x  _  e~ (°+j (n-n')P ^ 


(A.  6) 


Finally,  since 


1  _  e-(o+. j  (n+n' )) 
x  _  e~(a+j (n-n’)) 


2 


a  2  2 

(e  -  cosCn+n1))  +  sin  (n+iV) 

(e°  -  cos  (r|— n 1 )  )2  4-  sin2  (n-n1) 


e2  -  2  e  cosCri+n1)  +  1 

e2  -  2  e  cos(n-n’)  +  1 

cosh  a  -  cos  (r]+nf) 

cosh  a  -  cos  (rrn  * ) 


is  given  by 


^  1  i  .cosh  a 

G-  =  7“  log  ( - : - 

1  4  6  cosh  a 


cos (n+n1 )% 


-  cos (n-n’ ) 


As  G  tends  to  zero,  G^  becomes 


_  1  ,1  -  cos (n+n1), 

G1  ’  4  108 


(A.  8) 


(A. 9) 


For  all  (H1,?1)  in  a  small  neighborhood  of  (n,£),  the 
following  approximation  is  valid: 


82 


Re [log  (1  -  e 


(U-c 


+  j (n-n'))j j 


-  Re [log  (Ic-C’l  +  j(n-n'))] 


=  log 


J(j\-r]')2  +  (5-5' )2)  . 


(A. 10) 


Thus,  exhibits  a  logarithmic  singularity  of  the  form  (A. 10)  at 

(n',5')  =  (n,5). 
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Appendix  B 

Consider  the  series 


oo 

G  =  l  a  sin(nri)  sin(nn') 

^  _ n  n 

n=2 

A|5‘5’l  -„iwl 

e  e  1  1 

a  =  - - 

n  (3  n 

n 


6 

n 


Clearly 


I 

n=2 


1  <:  0  <  2  . 

¥ 


Put 


e-3(t)U-C'|  e-t|C-C 

a(t)  = - B(t) - t - 


3(t) 


■/ 


2  2 

t  -  a 


t  >  2 


\  -3(t)z  -tz 

F  (z ,  t)  =  e  -  e 


Then 


a(t) 


F(z,t)  dz  . 

\K-V\ 


(B.l) 


(B.2) 


(B.3) 


(B.4) 


(B.5) 


Since  t  <  $(t)  for  all  t  _>  2,  then  a(t)  >  0  for  all  t  2,  and 
certainly  so  is  n=2,3,...  . 
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That  {a  }  is  a  monotonically  decreasing  sequence  follows 


from  its  positiveness.  Since,  F  and  7^-  F  are  continuous  for  t  >  2 

3t  “ 

and  for  all  z  _>  |£-£T|,  and,  furthermore,  the  integrals 

»°°  rCJCt 

F(z,t)  dz  and  |  F(z,t)  dz  converge  uniformly,  then 

\Z-V\  |S-€*I 

[21,  Section  7-5] 


dt 


a(t)  = 


3 1 


F (z , t )  dz 


\z-v 


(  .« 


J 

ie-ef 


dz  ( 


F  (u, t )  du) 


(B.  6) 


with  the  help  of  (B.3)  and  (B.4).  Thus,  —  a(t)  <  0  on  [2,  00 ) . 

dt 

Consequently,  a(t)  is  a  monotonically  decreasing  function  for  all 
t  2,  and  certainly  so  is  a  ,  n=2,3,...  . 

00 

Thus,  G0  is  dominated  by  J  a  ,  a  series  of  positive 


n=2 


n 


monotonically  decreasing  terms.  Since 


6n+l  >  n  • 


n  =  2,3,  ... 


(B.7) 


then 


an  <  6  6  5  n=2»3» 


n 


n+1 


(B.8) 


Consequently 

N 


l  a  < 

n=2  n 


-32K-rl  -bn+1U-c 

e  e 


6, 


B. 


N+l 


(B.9) 


whence 
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I  S'  -  \  an  ^ 

n=2 


-32U-C 


(B. 10) 


-e, 


N+l 


since 


N+l 


goes  to  zero  as  N  goes  to  infinity. 


of  convergence  of  G^  is  therefore  exponential. 
Put 

b„-f  -i  ,  n-2,3,...  . 

n 


The  rate 


(B. 11) 


Then 


a  <  b  ,  n  =  2,3,...  .  (B.12) 

n  “  n 


Consequently 


00  oo 


(B. 13) 


which  can  be  proven  using  a  similar  procedure.  does  therefore  con¬ 

verge  uniformly  for  all  (r),£)  and  all  (r}f>£f)* 

All  the  above  results  can  also  be  deduced  from  Figure  B.l* 


Figure  B.l.  Pictorial  illustration  of  the  procedure  in 
the  Appendix. 
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Chapter  4 


MULTIPLE  CAPACITIVE  POSTS  IN  A 


RECTANGULAR  WAVEGUIDE 


Consider  a  system  of  posts  P  ,  P 


)  1  1  *  > 


located  close  to 


each  other  in  a  rectangular  waveguide.  These  posts  are  assumed  per¬ 
fectly  conducting,  of  arbitrary  shape  and  thickness,  and  uniform 
along  the  broad  side  of  the  waveguide,  i.e.,  of  the  capacitive  type. 
The  medium  filling  the  waveguide  is  assumed  linear,  homogeneous, 
isotropic,  and  dissipation  free,  and  is  therefore  characterized  by 
the  real  scalar  permittivity  £  and  the  real  scalar  permeability  y. 
The  problem  considered  is  depicted  in  Figure  1. 

1 .  Preliminary  Considerations 


Let  a  to  z  mode  of  unit  amplitude  be  incident  on  the 

posts  from  the  left.  This  mode  has  the  field  distribution 


> 


(1) 


,71  \ 

cos  (—  x)  e 
a 


where 
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z  =  0 


Figure  1.  p  capacitive  posts  in  a  rectangular  waveguide. 
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In  (2),  K  is  the  wave  number  of  the  waveguide  medium,  and  A  is  its 
wave  length.  Furthermore,  it  is  assumed  that  a  <  A  <  2a  and  2b  <?  A 
so  that  only  the  dominant  mode  can  propagate  in  the  waveguide. 

Since  each  post  is  uniform  along  the  x-axis,  and  since  the 
exciting  mode  has  no  x-component  of  electric  field,  neither  doe$  the 
scattered  field.  That  is,  the  scattered  field  is  TE  to  x,  and  can 
therefore  be  derived  from  an  electric  vector  potential  F  having  only 
an  x-coraponent  ip  [2,  Section  8-7]; 


F  -  i p*  . 


(3) 


The  scattered  field  is  given  in  terms  of  ip  by 


E  ~  -  V  x  ip  x 


HS  =  -7—  V  x  V  x  4>  x 

~  JCd\l 


(4) 


while  ip  itself  satisfies 


(V2  +  K2)lp  =  0  . 


(5) 


Expanding  (4)  in  rectangular  coordinates,  the  components  of  the 


scattered  field  are  found  to  be 
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(6) 


Furthermore,  since  each  post  is  uniform  along  the  x-axis, 

and  since  the  exciting  mode  has  an  x-component  of  magnetic  field 

7T 

that  varies  as  sin  (—  x) ,  so  does  the  scattered  field.  It  then 

a 

7T 

follow^  from  (6)  that  ip  must  contain  sin  (—  x)  as  its  x^dependent 

a 

factor.  Thus 


7T 

ip  =  sin  (—  x)  ^(y,z)  . 
a 

Substituting  (7)  into  (5)  and  psing  (2),  it  becomes 


(7) 


2  +  -^"2  ”  Yq)  >p(y»z)  =  o  .  (8) 

3y  9z 


The  components  of  the  scattered  field  are  now  given  by 
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ES  -  0 
x 

^S  /TT  N  3  I  /  x 

E  =  -  sxn  (—  x)  —  iHy,z) 
y  a  dz 

E  =  sxn  (—  x)  t-  ^(y,z) 
z  a  3y 

2 

s  ”V0  7T 

H  =  - -  sin  (—  x)  ip(y,z) 

x  joiy  a 


HS  = 

y  jwya 


cos  (—  x) 
a 


3y 


^(y,z) 


TT 

- -  COS 
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,TT  \ 
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_3_ 

9z 


^Cy^z) 


The  total  field,  incident  plus  scattered,  must  have  zero 
tangential  electric  field  at  the  waveguide  walls.  The  incident 
field  is  a  free  waveguide  mode,  and  does  therefore  have  zero  elec¬ 
tric  field  tangent  to  the  walls.  The  scattered  field  must  then  have 
zero  tangential  electric  field  at  the  walls.  This  is  readily  accom¬ 
plished  by  setting 


ip(y,z)  -  0,  y  =  0,b,  and  all  x  and  z  (10) 

dy 

which  are  the  required  boundary  conditions  on  the  scattered  field. 

The  boundary  conditions  (10),  once  satisfied  for  any  value 
of  x,  are  clearly  satisfied  for  all  values  of  x.  Thus,  the  problem 
is  basically  a  two-dimensional  scalar  one  that  can  entirely  be  worked 
out  in  some  x=constant  plane  within  the  waveguide.  Considerable 
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simplification  in  the  solution  can  result  from  choosing  the  x  *  ~ 
plane.  In  this  plane,  the  only  components  of  the  scattered  field 
are 


e5  ■  -  * 


s  ^0 

H  =  - -  \p(y9z) 

x  ju)y  r  w  ’ 


(11) 


whereas  those  of  the  incident  field  are 


E  =  e 


~V 


H  * 


~Y0  ~Yo2 


x  jojy 


(12) 


Finally,  an  important  fact  can  be  established  by  realizing 
that  the  set  of  equations  (8),  (10),  and  (11)  is  the  same  as  the  set 
satisfied  by  the  fields  in  a  parallel  plate  transmission  line,  ex¬ 
cept  for  Yq  replacing  jK.  That  is,  in  order  to  solve  any  problem 
of  capacitive  posts  in  a  rectangular  waveguide,  one  need  solve  only 
the  parallel  plate  transmission  line  problem  which  ha?  the  same 
cross  section,  at  the  same  time  replacing  jK  by  Yq*  In  the  latter 
problem,  the  exciting  field  is  a  TE^  to  x  mode,  again  with  Y„  re- 
placing  jK, 

No  further  simplification  is  possible,  but  i p  is  still  to  be 
determined.  This  can  be  accomplished  with  the  help  of  the  Green* s 
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function  for  TE 


In 


the  next  section. 


to  x  modes  in  a  rectangular  waveguide  obtained  in 


2. 


The  Green’s  Function  for  TE 
Waveguide 


In 


to  x  Modes  in  a  Rectangular 


Consider  a  magnetic  current  filament  M  directed  across  the 
waveguide  parallel  to  the  x-axis  and  located  at  (y*,  zf).  Further- 

7T 

more,  M  is  assumed  to  vary  as  sin  (—  x)  along  the  x-axis.  Figure  2 

a 

shows  the  situation  considered. 

Since  M  is  directed  along  the  x-axis,  the  field  produced 

must  have  an  x-component  of  magnetic  field  and  no  x-component  of 

7T 

electric  field.  Furthermore,  since  M  varies  as  sin  (—  x)  along  the 

a 

x-axis,  so  does  H  .  This  can  readily  be  verified,  for  instance, 
x 

from  the  reciprocity  theorem  [2,  Section  3-8].  An  electric  vector 

potential  having  only  an  x-component  proportional  to  H  ,  as  is 

x 

seen  in  Section  1,  can  then  be  used  to  derive  all  field  components. 

Only  TE^  to  x  modes  can  be  excited  in  the  waveguide, 

since  these  are  the  only  modes  having  an  H  component  that  varies 

x 

7T 

as  sin  (—  x)  and  no  E  component  [2,  Section  4-4].  The  potential 
a  x 

function  ip  due  to  the  filament,  relabeled  G,  is  therefore  referred 
to  as  the  Green’s  function  for  TE  to  x  mo^es  in  a  rectangular 
waveguide.  Below,  G  is  found  as  a  series  of  these  modes. 

Put 


G  =  sin  (—  x)  G(y , z) 
a 


(13) 


G(y,z)  then  satisfies,  for  each  x»  the  wave  equation 
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Figure  2.  A  magnetic  current  filament  M  in  a  rectangular 
waveguide  parallel  to  the  x-axis. 
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2  2 

C— 2  +  “^“2  "  G^y,z^  =  ~  5 (y-y * >  6(z-z')  (i4) 

3y  3z 


together  with  the  boundary  conditions 


3y 


G(y,z)  =  0 


y  =  0,b,  and  all  z. 


(15) 


In  (14),  6  is  the  Dirac  delta  function,  Myltiplying  both  sides  of 


,n7T 


(14)  by  cos  (—  y)  ,  then  integrating  over  y  from  0  to  b,  it  becomes 


,  d  2  /HTT.2 

(T"2  ~  Y0  ~  <"b  >  ) 

dz 


.nir 


G(y,z)  cos  (—  y)  dy 

D 


,nTT 


cos  (—  y ? )  6  (z-z 1 ) 
b 


(16) 


Put 


Gn(z)  = 


G(y,z)  cos  y)  dy 


,  .n7i\  2  2 

)  +  Yr 


,  n  =  0,1,2, 


(17) 


The  one-dimensional  wave  equation  (16)  then  becomes 


<72  -  G„<’>  ■ 

dz 


-  cos  (”^-  y')  6(z-zf) 


(18) 


For  the  solution  of  (18)  to  represent  waves  traveling  away 


from  the  source,  G  must  be  of  the  form 
n 


Gn(z)  -  < 


-Y  z 

A  e  n 
n 

z  >  z 1 

Yz 
_  n 

B  e 

z r  >  z 

(19) 


where  A  and  B  are  constants  to  be  determined.  Since  G  is  pro- 
n  n 

portional  to  H^,  it  is  continuous  across  the  filament  at  ?  =  z% 

[2,  Section  1-14],  and  so  is  G^.  Thus 

-y  z1  y  z* 

A  e  n  -  B  e  n  =0.  (20) 

n  n 


Furthermore,  integrating  (18)  over  z  from  z?  -  A  to  z1  +  A,  then 
letting  A  go  to  zero,  there  then  results 


dz 


G  (z) 
n 


/HTT  ?  \ 
=  -  cos  (-T-  y  ) 


(21) 


That  is, 
Thus 


—  G  is  discontinuous  at 
dz  n 


z 


z 1  by  the  amount 


/UTT  | 
-cos  (— y 


-Y  z 1  y  z 1 

n  n 

A  Y  e  +  B  y  e  ^cos 

n  n  n  n 


(22) 


Solving  (20)  and  (22)  simultaneously,  A^  and  are  found 


to  be 


*  -L-  ,UM  t  x 

A  =  “5"  cos  (—  y1)  e 
n  2y  b 

'n 


Y  z 1 

*  n 


(23) 


i  ~Y  z ' 

„  1  /HIT  ,  ^  n 

Bn  '  W  C°S  CT  y  )  6 
n 


(24) 


Combining  (23)  and  (24)  with  (19),  becomes 


Gn(z)  =  2TCOS  (Ty,)  e 
'n 


-Y  |z-zf 


,  n=0, 1, 2, ...  .  (25) 


By  Fourier  theory  [15,  Section  43],  (17)  can  be  inverted  as 


where 


G(y , z | y 1 ,zf)  = 


t  Z‘. 

n=0 


„  — Y  z-z ' 

-nTT  .  .nir  . .  'n' 

cos  (—  y)  cos  (-7-  y')  e 
b  b 


2Y, 


n 


(26) 


£ 

n 


1 

2 


n  «  0 

n  >  1 


Clearly,  G  satisfies  the  boundary  conditions  (15). 


(27) 


3.  Basic  Formulation 

As  is  pointed  out  in  Section  1,  the  problem  is  a  two- 

a 

dimensional  scalar  one  that  can  entirely  be  worked  out  in  t\ie  x  =  ~ 
plane  within  the  waveguide.  Thus,  all  source  and  field  points  are, 
hereafter,  assumed  located  in  this  plane.  The  incident  field  is 
then  given  by  (12),  whereas  that  scattered  from  the  posts  is  given 
by  (11). 

Let  (E^,  jH^)  be  incident  while  all  the  posts  are  absent,  and 

(ECJ),  H(J))  be  the  field  produced  by  an  electric  current  of  density 

®  Jm9  where  Jm  is  the  current  on  C™ ,  while  all  the  posts  are 
m=l  i 

absent.  By  the  uniqueness  theorem  [2,  Section  3-3],  (E  +  E(J) , 

H1  +  H(J))  is  identical  with  the  original  field  whenever 

nm  x  (E1  +  E(J))  =0  on  Cm  .  (28) 

In  (28),  nm  is  the  outward  unit  vector  normal  to  Cm  and  tangent  to 
the  x  =  -z  plane.  (JiQ)>  HQ.))  must  then  have  the  field  distribution 


(11) .  Since 
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m 


n  x  (H(J) 


-  H(J) 


v=0 . 


)  = 


Tm  m 

J  on  C 


(29) 


v=0 


where  v  is  the  distance  along  n™  from  Cm,  and  H(J)  has  only  an  x- 
component  that  does  not  vary  with  x,  _Jm  has  only  an  x-transverse 
component  Jm  that  does  not  vary  with  x: 


t111  Tm/  \  *-ra  <_m 

J  =  J  (y,  z)  t_  «  J  (t)  t 


(30) 


In  (30),  _tm  is  the  counterclockwise  unit  vector  along  Cm,  and  t  is 
the  distance  along  Cm  from  an  arbitrary,  but  fixed,  point  to  (y,z). 
Consider  the  magnetic  current  distribution 


V’  x  Jm(t’)  6(v')  tm  on  Cm  .  (31) 

Y0 


Since 


f— » *  i  m  r—i  *  1  m 

V'  x  (J)  _t  =  V ' X  t 


3v»  v  - 


cj)  n®  x 


in 


3v*  - 


(32) 


for  any  function  (})  independent  of  x,  (31)  becomes 


-  2 
Y0 


Jm(t') 


dVT 


6(vf) 


on  C 


m 


(33) 


The  collection  ^  produces  a  field  identical  with  (E(J) »  H(.J)) 
m=l 

[22].  By  definition,  G(y,z|y,,zl)  x  is  the  electric  vector  poten¬ 
tial  produced  at  any  point  (y,z)  by  a  unit  magnetic  current  fila¬ 
ment  in  the  x-direction  located  at  (y1,  zf).  Thus,  by  superposition, 
produces  at  (y,z)  the  electric  vector  potential  ^ x ,  where 
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! 

v2 


<K  y»z) 


=  JMi 


Y0  w-1 


Ay',*’)  [ 


m 


-v’ 


/ 


dt 1  =  / (dy ' )  +  (dz ' ) 


i  d  ^ 

G(y,z|y',z')  6 (v ’ )  dv']dt' 


(34) 

J 


In  (34),  G  is  given  by  (26)  and  (27),  and  primed  and  unprimed  coordi¬ 
nates  denote,  respextively ,  source  and  field  points.  Whence 


41  =  r-M!  I 

r 0  m=;L 


m 


jI"  G(y,z|y',z')  dt'  . 


(35) 


Since  G(y,z|y!,z!)  is  a  solution  of  the  homogeneous  wave  equa¬ 
tion  (8)  for  all  (y,z)  ^  (y^z1),  so  is  \\).  Furthermore,  ^  satisfies 
the  boundary  conditions  (10)  by  virtue  of  (15).  Thus,  \p9  and  con¬ 
sequently  the  complete  field  solution  of  the  problem,  can  be  found 
once  all  Jm  are  known.  Using  (4),  (12),  and  (35),  (28)  becomes 


n*  x  [e‘V 


Y0  m*1  C 


ra 


Jm(y',zf)  G(y,z|y'  ,z')  dt')x] 


=0,  (y,z)  e  C  ,  1  <  £  <  p  . 


(36) 


Since 


£  .  .  £  ,  3  .  £ 

_n  x  (V  x  (p  x)  =  n  x  (V  <p  x  x)  =  n  x  (—  <j>  n  x  x) 


3v 


<j>  x 


(37) 


(36)  can  be  put  in  the  form 
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sin(0 (n  ,2))  e 


'V  .  jMd 


Yn  rn  1 


Jm(y' ,z*) 


,ra 


G(y > z | y *  # z 1 )  dt’ 


*  0,  (y,z)  £  c\  1  £  &  p  (38) 

l 

which  is  an  integral  equation  for  J.  In  (38),  0(n  is  the  angle 

Z 

n  makes  with  the  y-axis  at  (y,z). 

The  higher  order  (n  _>  1)  modes  excited  are  evanescent,  i.e,, 
decay  exponentially  with  distance  from  the  posts.  Thus,  at  distances 
sufficiently  far  from  the  posts,  only  the  dominant  mode  (n=0)  can 
exist.  The  reflection  coefficient  of  the  dominant  mode  is  readily 
found  from  (11),  (12),  (26),  (27),  and  (35)  as 


r  =  2m  V 

'  2V  m=l 


Jm(y',z')  sin  (6(nm,£))  e 


"V 


dt'  . 


(39) 


m 


The  transmission  coefficient  of  the  dominant  mode  is  then 


T  =  1  + 


2V 


m=l 


m 


sin(0(n  ,^'))  e 


V 


dt' 


(40) 


m 


4.  The  Scattering  Matrix 


Following  Montgomery  et  al .  [10,  Section  5-14],  the  scatter¬ 
ing  matrix  of  the  posts  is  defined  as 


S  = 


(41) 
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In  (41),  and  are,  respectively,  the  amplitudes  of  the  domi¬ 

nant  mode  reflected  to  the  left  and  transmitted  to  the  right  of  the 
posts  due  to  an  incident  TE^q  to  z  mode  of  unit  amplitude  from  the 
left.  Consequently,  and  are  given  by  (39)  and  (40),  respec¬ 
tively  . 

Similarly,  and  are’  resPectively ,  the  reflection  and 

transmission  coefficients  of  a  TE^q  to  z  mode  of  unit  amplitude  inci- 
dent  on  the  posts  from  the  right.  In  the  x  -  ~  plane,  this  mode  has 
the  field  distribution 

(42) 


i  V 

E  =  e 

y 


x  jwp 


V 


The  previous  analysis  carries  through  in  this  case, 
the  scattered  field  is  given  by  (11)  and  (35) ,  but  with  ^J1 


Thus, 


ra=l 


now  replacing  in  (35),  determined  by  solving  the  integral 
equation 


sin  (6 (n^,y) )  e 


■v 


_ 

2 


P 

l 


Yn  m=1 


J,m(y\z') 


m 


9v3v* 


0 (y , z I y 1 , z ' )  dt' 


=  0  ,  (y,z)  e  C 1  £  £  £  p  (43) 

rather  than  (38).  It  then  follows  from  (11),  (26),  (27)  and  (35) 


that 


102 


22 


=  Jhl 


2Y0b  m=l 


Y  z 

J'ra(y'  ,z')  sin  (6(rim,x'))  e  ^  dt 1  (44) 


m 


s  =  !  +  Ml.  y 

°12  x  9v  h  L 


2y0b  m=l  ; 


— Y  z  1 

,z')  sin  e  ^  dt' 


(45) 


,m 


The  scattering  matrix  is  both  symmetric  and  unitary.  That  is 


S  =  S 


(46) 


ssH  =  sHs  =  u 


where  T  and  H  denote  matrix  transpose  and  Hermitian,  respectively, 

and  U  is  the  identity  matrix. 

11  2  2 

Let  (E  ,  H  )  and  (E  >  H  )  be  the  z-transverse  fields  in  t\\e 
waveguide,  sufficiently  far  from  the  posts,  due  to  t0  z  m°des 

of  arbitrary  amplitudes  c^  and  c^  incident  from  the  left  and  from 
the  right  of  the  posts,  respectively.  It  then  follows  from  (11), 


(12),  (26),  (27),  (35),  (39),  (40),  (42),  (44),  and  (45)  that 


e1  =  <; 


“V 


Ynz 


C1  (e  +  sn  e  >2 


q  “V 

C1  S21  6  * 


z  «  0 


z  »  0 


H1  = 


"Ynz  Yoz 

-n  c  (e  -  S  e  )  x  2  «  0 


_n0  cl  s21  6  2 


z  »  0 


(47) 
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C  rj  S.  «  G 

2  12  2. 

Y0Z  'V 
c2  <e  +  S22  6  }  * 

V 

n0  C2  S12  6  - 

,  V  e  'V, 

nO  C2  **  S22  6  ^  - 


z  «  0 


z  >>  0 


z  «  0 


z  »  0  . 


In  (.47),  rig  is  the  characteristic  admittance  of  the  dominant  wave¬ 
guide  mode: 


.  j_  =  A 
no  ”  c0  ~  jwy 


(48) 


Let  W  be  the  curve  enclosing  the  surface  area  of  the  waveguide 
in  the  x  =  —  plane  between  the  z  =  and  lines,  for  some  <<  0 
and  »  0.  The  reciprocity  theorem  then  states  that 


(E1  x  H2  -  E2  x  H1)  •  n  dt  =  0  (49) 

W 

where  n  is  the  outward  unit  vector  normal  to  W  and  tangent  to  the 
x  =  ^  plane.  Substituting  (47)  into  (49),  there  then  results 

2b  cx  S12  c2  =  2b  ^  Sn  c2  (50) 


whence 


S12  ”  S21 


(51) 


The  scattering  matrix  is  symmetric  if  and  only  if  the  whole  structure 
is  reciprocal. 
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That  S  is  unitary  follows  from  conservation  of  power  con¬ 
siderations.  Let  the  two  dominant  modes  be  simultaneously  incident 
on  the  posts  from  the  left  and  from  the  right.  The  complex  power 

a 

scattered  far  to  the  left  and  to  the  right  of  the  posts  in  the  x  *  ~ 
plane  is  basically 


P 

sc 


bnr 


(IC1S11+C2S12I 


lClS21  +  C2S22' 


(52) 


l 


whereas  that  incident  is  given  by 


'in-^O  <lcll2+  M2' 


(53) 


Since  the  structure  is  lossless,  and  since  P.  and  P  are  real,  they 

m  sc 

must  be  equal.  Put 


(54) 


Then 


+H+  -^H  H 

bnQ  c  c  =  br|Q  c  S  S  c 


(55) 


or 


sHs  =  u  . 


5.  The  Impedance  Matrix 

Let  TE^q  to  z  modes  of  arbitrary  amplitudes  c^  and  c^  be 
incident  on  the  posts  from  the  left  and  from  the  right,  respectively. 
Let  and  v^  be,  respectively,  the  amplitudes  of  the 
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component  far  to  the  left  and  to  the  right  of  the  posts  referred 
to  the  z=0  plane.  It  then  follows  from  (47)  that 


V1  (1  +  Sll)  C1  +  S12  C2 


(57) 


V2  =  S21  C1  +  (1  +  S22)  C2  ‘ 


(58) 


The  choice  of  z  =  0  as  a  reference  plane  is  only  a  matter  of  con* 
venience .  In  matrix  form  (57)  and  (58)  become 


v  =  (U  +  S)  c 


(59) 


where 


v  = 


2 

L.  J 


(60) 


Similarly,  let  and  i^  be,  respectively,  the  amplitudes  of 
the  component  far  to  the  left  and  to  the  right  of  the  posts 
extrapolated  back  to  the  z  =  0  plane.  Then 


C0  11  (1  Sll)  C1  S12  C2 


(61) 


C0  12  ~  S21  C1  +  (1  S22^  C2 


(62) 


In  matrix  form,  (61)  and  (62)  become 


i  =  (U  -  S)  c 


(63) 


where 


1C6 


To  relate  to  network  theory,  let  (v^,  -1^)  and  ( v 12)  be 
the  complex  voltage-current  pairs  at  the  terminals  of  a  two-port 
network  {10,  Section  4-5].  Then 


v  =  Z  i  (65) 

where  Z  is  the  network  impedance  matrix.  From  (59)  and  (63),  Z  is 
readily  found  as 

Z  =  (U  +  S)(U  -  S)"1  (66) 

Since  S  is  symmetric,  so  is  Z.  Furthermore,  since 


Z  =  £0(U  +  S)(U  -  S)  1  =  c0(sHs  +  S)(SHS  -  S)  1 

=  c0(sH  +  U)(SH  -  U)"1 
=  -  ZH  (67) 


the  elements  of  Z  are  pure  imaginary.  Thus 


Z  =  j 


X 


11  12 


21  22 


(68) 


X12  "  X21 
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6 .  The  Equivalent  Network 


The  complete  field  is  seldom  needed.  Rather,  the  effect  of 

the  posts  on  the  dominant  waveguide  mode  is  what  must  accurately  be 

discribed.  From  an  engineering  perspective,  a  description  in  terms 

of  a  network  of  lumped  elements  is  preferred. 

The  effect  of  the  posts  on  the  dominant  waveguide  mode  is 

fully  described  by  the  impedance  matrix  Z.  Such  a  representation  can 

be  realized  in  the  form  of  a  two-port  T-network  [10,  Section  4-5]. 

The  characteristic  impedances  of  the  TE.,  to  x  modes  are 

In 

given  by 


«n  '  - 


'(•El)2  +  #): 2 

b  a 


(1,2 

a 


n  >  1 


c  - 


> 


(69) 


Since  these  modes  are  evanescent,  the  energy  stored  close  to  the 
posts,  in  view  of  (69),  is  predominantly  electric.  This  effect  can 
suitably  be  represented  by  a  capacitor  in  the  shunt  arm  of  the  net¬ 
work.  The  elements  in  the  series  arms  are  also  capacitors  to  account 
for  the  change  difference  across  the  posts  in  the  z  direction.  The 
equivalent  network  of  the  posts  is  shown  in  Figure  3. 


7 .  Solution  of  the  Intergal  Equation 


The  integral  equation  (38)  can  be  written  in  compact  form  as 
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j(X|i  x12) 
— )\— 


i 


o  i 


z  =  0 


j  (  Xg2  X  |2 ) 

- 


jx 


12 


z  =  0 


Figure  3.  The  equivalent  network  of  the  posts. 
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l  Zm  (Jm)  =  sin  (0(nL  £))  V 
m=l 


zm  (Jm}  =  ^ 

Yo  cm 


jm(y'5z')  'dvdv7  G^y>zly'’z')  dt' 


~Y0Z  £ 

V  =  e  ,  (y ,  z)  e  C  ,  1  £  £  £  p 


(70) 


An  exact  solution  of  (70)  can  rarely  be  obtained  and  an  approximate 
solution  has  then  to  be  sought. 

Let  each  C™  be  approximated  by  a  polygon  =  {s™,  ,  •  .  *  ,  Smm  } 

as  shown  in  Figure  4,  and  put 


Jm  (y\zf) 


m 

j=i 


ro  Tm  j  y  i  \ 

i.  Jj(y  »z  )  • 


(71) 


In  (71),  I™  are  complex  coefficients  to  be  determined,  whereas  each 

Jm  is  a  real  function  that  vanishes  on  all  ^  ra,  but  is  otherwise 
J 

unspecified.  Thus,  0  (nt  ,  y_)  is  constant  on  S™: 


0(nm,  z)  =  0(nm,  Z)  =  0m  on  S™  . 

—  '**■'  — j  j  j 


(72) 


Substituting  (71)  and  (72)  into  (70) ,  there  then  results 


l  f  i! 


m=l  j=l 


m  m  t  -.Ills  .  .  /  TT 

Z  (J  )  +  r  =  sm  (0^)  V  , 


( y » z )  £  S^,  l£££p,  l£i<  q^.  (73) 


The  integrals  in  (73)  are  taken  over  S^,  and  r  is  a  residual  term. 
A  Galerkin  solution  [4,  Section  1-3]  can  be  obtained  by  requiring 


uo 


) 


Figure  4. 


m  v 

C  approximated  by  a  polygon  l 


m 


Ill 


m 


that  r  be  orthogonal  to  all 

Define  the  inner  product 


<A,  B>  = 


* 

A  B  dt 


(74) 


£=1 


where  *  denotes  complex  conjugate.  Taking  the  inner  product  of 

l 

(73)  with  each  JL  and  enforcing  the  Galerkin  condition 


<r ,  =  0 

x 


1  <  SL  <  P  *  1  _<  i  £  q 


(75) 


it  becomes 


m 

v 


)  <im  zm  (jm),  j£> 
m=l  j=l  J  J  1  1 


o  i 

sin  (0p  <V,  J*>  , 


1  <  £  <  p  ,  1  <  i  <  q 


(76) 


The  system  of  equations  (76)  can  be  put  in  the  matrix  form 


Z  I  =  D  V  . 


(77) 


In  (77) ,  Z  is  a  p  by  p  block  matrix  whose  £mth  block  is  the  matrix 


z*"  -  [z*m]  .  -  [<z"  u”),  j'->]  , 


(78) 


D  is  a  p  by  p  block  diagonal  matrix  whose  ££th  block  is  the  diagonal 
matrix 


_££  ,■_££-,  r  .  ,q£m 

D  -  ‘“id  i  z-‘sm  (0i)!  i  l 

q  xq  q  xq 


(79) 


and  I  and  V  are  p  segment  vectors  whose  mth  and  £th  segments  are 


the  vectors 
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T1  =  Iimj 

J  qrax! 


(80) 


V 


[vjj  -  [<V,  A]  £ 

q  X1  q  xl 


(81) 


respectively. 

The  currents  J™  given  by  (71),  with  the  coefficients  I™ 
determined  from  (77),  form  the  Galerkin  solution  of  (70).  A  galerkin 
solution  of  (43)  can  be  obtained  in  a  similar  manner.  Clearly,  then, 
using  the  same  J™,  the  solution  is  given  by  (71),  but  with  coefficients 
now  determined  by  solving  (77)  with  the  right-hand  side  vector  V  con¬ 
jugated 


8 .  Evaluation  of  the  System  of  Equations 

The  construction  of  Z  in  (77)  constitutes  a  large  portion  of 
the  work  involved  in  the  numerical  solution.  An  efficient  evaluation 
of  the  elements  of  Z  is  therefore  necessary  for  its  success. 

The  ijth  element  of  the  &mth  block  is  given  by 


7Vm  j  coy 
2 


Yi 


Ji  (y , z)  dt 


Jj  3~v~a~vr  G(y>z\y' >z")  dt' 


o  s: 

1 


sm 

3 


(82) 


m 


where  are  so  far  unspecified.  A  particularly  simple  choice  for 

Tm  . 

J .  is 


X13 


Jj(y'z’)  =  ' 


m 


1  (y '  »z')  e  S 


o  (y'  »z')  C  S*  “  ® 


(83) 


which  corresponds  to  a  pulse  expansion  of  J™.  Z^J1  then  becomes 


,£m  _  juy 
ij  '  v2 


dt 
J  9,  i 

Y0  S.  S. 

l  J 


a2 


m 


9v3v 


7  G(y , z |y ' , z  ' )  dt'  . 


(84) 


Put 


Ry(y.z)  = 
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'  m 


9v9v 


T  G(y,z|y'  ,z')  dt',  (y,z)  e  S.' 


£ 


(85) 


By  the  first  mean  value  theorem  of  integration  [16,  Section  7-18], 

£ 

there  exists  a  point  (yo»zQ)  e  such  that 


„£m  jooy  £  £m  ,  , 

Z.  .  =  L.  R.  .  (y  ,z  ) 
ij  2  l  ij  o  o 

Y0 


(86) 


where 


T  £  f,  £ 

h  J 


H.2  .  ,  H 

yi>  +  <zi+l 


£.2 

zi} 


(87) 


is  the  length  of  S^. 
Put 


G (y , z | y ' , z ' )  =  (Gq  +  +  g2) (y,z|y' ,z')  . 


(88) 


In  (88) 
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G0(y,z|y*>Z') 


-1  blz"z  i^o 

j2iC 


nTT  | 


G1(y»z|y'  .z')  =  \  l  ^  COS  (^—  y)  cos  y')e  b 

n=l 


z-z 


G2(y,z|y' ,z') 


uu 

1  r  ytlTT  v  /n7T  ,  . 

-  )  cos  (—  y)  cos  (—  y’) 

TT  b  b 

n-1 


TT  i 


nTT  i 


( 


b  z-z '  b 

e  e 


z-z 


n 


(89) 


Y, 


0 


•  71  o 

J  b  1  0 


-nir.  2  ,  2 

<T>  +Yo 


TT 

b 


3  . 

n 


(90) 


The  decomposition  (88)  therefore  amounts  to  expressing  the  dynamic 
Green’s  function  G  in  terms  of  a  dominant  mode  wave  G^,  the  corre¬ 
sponding  static  Green’s  function  G^,  which  can  be  obtained  from  (26) 
by  dropping  the  n  =  0  term  and  setting  equal  to  zero  in  the  re¬ 
maining  terms,  plus  a  correction  series  G^. 

The  series  defining  G^  is  readily  summed  to  give 

G1(y,z|y',z')  =  ±\z- z' | 

-  log [ cosh I  z-z ' |)  -  cos  (y+y'))) 


x  (cosh  (-^  |  z-z'  | )  -  cos(-g  (y-y')))]  (91) 
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where  log  denotes  the  natural  logarithm.  The  details  of  the  summation 
are  given  in  Appendix  A.  The  series  in  G^  is  dominated  by  an  exponen¬ 
tially  convergent  series  of  positive  monotonically  decreasing  terms,  as 
is  shown  in  Appendix  B  of  Chapter  3,  and  can  therefore  be  summed  di¬ 
rectly  at  a  minimal  cost. 

The  normal  derivatives  of  G  at  (y  , z  )  e  S . ,  and  at  any 

o  o  i 

(y ' , z ' )  e  can  be  computed  numerically  through  finite  differences. 
Thus,  for  all  (£,i)  4  (m,j). 


M  M' 


9v3v 


7  G(yo,zo|y',z') 


cc 


l  J  ,Vw<VV<J2HV‘Jyk’-,k')' 

c=l  k  =1 


(92) 


Here,  (y  ,  z  )  are  the  pivot  points  along  the  normal  to  S  at  (y  ,z  ), 
k  k  i  o  o 

and  b,  are  the  coefficients  and  c  is  the  multiplier  factor  of  the  Mth 
k 

difference  formula.  The  parameters  for  the  M'th  difference  formula, 
(y^,,z^,),  b^,,  and  c',  are  similarly  defined. 

When  dealing  with  the  diagonal  elements  of  Z,((£,i)  =  (m,j)), 

G ^  offers  a  logarithmic  singularity  at  (yo>zQ)  that  requires  particular 
attention.  In  Appendix  A,  the  singular  part  of  is  found  to  be 


G1  (y  ,z  |y',z')  =  — y  log[.2(cosh(£  (z  -z')) 
Is  o  o 1  4tt  bo 


-  cos  (•£  (yo~y  ’ ) ) )  ]  .  (93) 


Put 


Glp(yo’Zo|y'’Z,)  =  (G1  ■  Gls)(yo’Zoly,>Z,)  • 


(94) 
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Then 


G(>'o-zoly'-'') 


0 

j26  : 

1  /TT\  „  r  e 

oT  Or)  Re[ - ~ - 

8b  b  ,  ,  2  ,u  ox 

sinh  (-  -) 


+  3^  <V<VG2)(VZoly'>z') 


1  (?)  net  e 


0 

j  20^ 


8b  b  .  ,  2  ,n  os 

sxnh  (-  -p 


] 


M  M' 


+  kl1bkbi'(G0«lp+<!2)(yk>2klyrzk') 


where 


(95) 


a  =  (zo-z’)  +  j  (yo-y')  • 


(96) 


The  truth  of  the  first  term  on  the  right-hand  side  of  (95)  can  readily 
be  verified  in  a  straightforward  manner  by  carrying  out  the  normal 


differentiations,  and  Re(z)  denotes  the  real  part  of  z. 

2 

The  evaluation  of  Z . .  is  now  completed  by  integrating  — r  G 

lj  dVoV 


m 


over  S_.  ,  and  for  that  purpose  any  quadrature  rule  can  be  used.  Thus, 


T  £  T  111  M 

n  •  L» .  L ,  N 

im  =  jwy  i  j  r 

ij  v  2  2  \ 

Y  Q  u-1 


\ 


M  M' 


2  J  ,WW’rtl>;k''4'>'  •  (97) 

k=l  k-1 


x 
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In  (97),  (y^,,z^,)  are  the  pivot  points  of  the  M'  difference  formula 

along  the  normal  to  S?  at  ((1-p  )y™  +  p  y? ,  ,  (1— p  )z?  +  p  z™ , 1  )  ,  and 

3  u  3  u  3+1  ru  3  ru  j+1 

q  are  the  coefficients  and  p  determine  the  location  of  the  abscissas 
u  u 

of  the  Nth  integration  rule.  The  diagonal  elements  of  Z  are  then  given 
by 


z££  _  jioy  l£  rgM 

ii  2  i  1  ii 


u=l  k=l  k  =1 

(98) 


where 


U  1_  .TT. 

Gii  8b  V 


o 

320. 

l 


Re  [- 


0  .  ,  2  ,U  O, 

S£  Sinh  (b  2} 
i 


]  dt’  . 


(99) 


The  singularity  in  the  integrand  at  (y  ,zq)  is  not  integrable.  However, 
the  limit  of  the  integral  as  (y,z)  approaches  (y  ,z  )  along  the  normal 
at  (y  ,z  ),  not  its  value  at  (y  , z  ),  is  what  is  actually  needed.  The 
integral  is  elementary,  nevertheless,  and  is  readily  evaluated  to  give 


£  £  £ 

G££  =  Re[e”  1  (coth(y  -  e  1)  +  coth  (y  y  e  1))]  (100) 

ii  4b  b  Z  o  l 


Z  Z 

where  L  is  the  distance  from  (y^  z to  (yQ>  ZQ)  • 


The  ith  element  of  the  Zth  segment  of  V  is  given  by 
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v£  = 

l 


i  "V 

J^(y » z)  e  dt 


which,  upon  using  (83),  becomes 


v£  = 

1 


"V 

e  u  dt  . 


(101) 


(102) 


The  Integration  in  (102)  can  be  carried  out  exactly.  However,  there 
exists  a  point  (yo,ZQ)  such  that 


T72  T2  Y0Zo 
V.  =  L.  e 


(103) 


1  i 

Actually,  finding  such  points  (y  ,  z  )  and  (y  ,ZQ)  is  at 

least  as  difficult  as  computing  the  integrals  themselves.  For  suf- 

£  2 

ficiently  small  L^,  however,  the  mid-point  of  can  replace  these 
points  while  introducing  only  very  little  error.  The  system  of  equa¬ 
tions  thus  obtained,  clearly,  is  one  that  results  from  enforcing  the 
point  matching  condition 


r (y , z)  =  0  , 


(y,z)  e  {( 


2^2 

yi+i  +  yi 


t  .  t 

z.,t  +  z.  0 

—  1  2 - 31  |l  £  2  £  p,  1  £  i  £  q  ) 


in  (73) ,  rather  than  the  Galerkin  condition  (75) . 


(104) 


9 .  Numerical  Results 

The  solution  procedure  presented  is  readily  translated  into  a 
computer  program.  The  elements  of  the  scattering  matrix  and  the 
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reactances  of  the  equivalent  T-  or  II-network  are  basically  the 
parameters  to  be  computed. 

The  scattering  parameters,  thanks  to  (71),  (72),  and  (83),  are 
computed  by 


 jwy 


m 

i  \ 


l  /  I.  sm  0, 

11  2^nb  J  rj 


"V 

e  U  dt' 


0  ra=l  j=l 


sm 

3 


S21  "  1  + 


m 

jwy  ^  \  .  /nmx 

27b  2  '•  Ii  sin(0i)  , 

'0  711=-- 1  j=l  J  1  ; 


V’  ,7. 

e  dt 


m 


m 

j wy  ^  \  Tfm  .  ,0mN 

i-  i.  sin(fV 

22  2V  IFlj=l  ]  3  J 


■Ynz 


sm 

3 


Sl2  =  1  + 


3<4J  ' 

2y„b  ^  . L,  ~j 
'0  m=l  3=1  J 


\  i:m  since?) 

L  1  3  J 


dt' 


"V 


dt' 


sra 

3 


(105) 


where  im  and  I !m  are  the  solutions  of  (77)  with  right-hand  side  vectors 
3  3 

V  and  V*,  respectively.  The  impedance  matrix  is  then  computed  through 
(66).  In  carrying  out  the  integrations  in  (105),  and  also  in  (97)  and 
(98),  an  eight-point  Gauss-Radau  quadrature  rule  [17]  is  used.  The 
parameters  of  the  integration  rule  are  shown  in  Table  1  of  Chapter  3, 
For  differentiation,  a  symmetric  finite  difference  formula  with 
M  =  2,  b^  =  1,  and  b^  =  -1  is  used.  The  multiplier  factor  c  is  set 
equal  to  two  tenths  of  the  length  of  the  segment  where  differentia¬ 
tion  occurs. 

Because  of  the  approximations  involved  in  the  solution,  the 
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scattering  matrix  need  no  longer  be  symmetric  nor  unitary.  To  deter¬ 
mine  the  impedance  matrix,  and  are  first  replaced  by  their 

average 


S 

av 


2  (S12  + 


S21) ' 


(106) 


The  impedance  matrix  can  then  have  a  non-zero  real  part. 

To  test  the  solution,  the  computer  program  is  run  for  a  few 
selected  problems.  In  particular,  the  problems  of  the  circular  post 
and  of  the  symmetrical  thin  strip  are  considered.  Some  of  the  results 
obtained  are  plotted  in  Figures  5-8, 

In  all  the  cases,  the  convergence  for  capacitive  susceptance 
B^  is  monotonic  and  from  above,  as  can  be  seen  for  the  thin  strip  in 
Figure  8.  The  computed  susceptances  are  found  to  agree  well  with  the 
data  in  the  Waveguide  Handbook  [3],  with  only  a  few  segments  needed 
even  for  large  posts.  A  complete  assessment  of  the  solution  should 
also  consider  the  (Frobenius)  norm  of  the  real  part  of  Y  =  Z  ^  and 
the  modulus  of  difference  in  transmission  coefficients.  These  two 
numbers  are  computed  in  all  program  runs,  and  are  usually  0(10  ^). 

10,  Concluding  Remarks 

The  system  of  capacitive  posts  in  a  rectangular  waveguide,  i.e,, 
of  metallic  obstacles  that  are  uniform  along  the  broad  side  of  the 
waveguide,  but  are  otherwise  of  arbitrary  shape  and  thickness,  has 
been  considered  in  this  chapter. 
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Figure  5.  The  series  susceptance  of  the  centered  circular  post. 
The  number  of  segments  used  is  30. 
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Figure  6 


The  parallel  susceptance  of  the  centered  circular  post, 
number  of  segments  used  is  30. 


The 
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B  X0 


Figure  7.  Network  susceptance  of  the  symmetrical  thin  strip.  The 
number  of  segments  used  is  20. 
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B  Xq 


NO.  OF  SEGMENTS 


The  convergence  behavior  for  the  symmetrical  thin  strip 

(f^  =  0.6,  £  =  0.5). 

Aq  b 


Figure  8. 
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A  complete  field  analysis  of  the  problem  is  first  given.  From 
the  analysis,  the  scattering  and  impedance  matrix  representations 
of  the  system  of  posts  fully  describing  its  effect  on  the  dominant 
waveguide  mode  are  obtained.  Since  the  whole  structure  is  both  recipro¬ 
cal  and  lossless,  the  scattering  and  impedance  matrices  are  symmetric 
and  unitary,  and  symmetric  and  pure  imaginary,  respectively.  The 
latter  is  then  realized  in  the  form  of  a  T-network  of  capacitive  ele¬ 
ments  • 

The  susceptances  of  some  post  configurations  are  computed  in 
Section  9.  In  the  actual  computation,  pulse  expansions  of  the  currents 
induced  on  the  posts  are  used.  Although  chosen  primarily  so  as  to  ren¬ 
der  the  procedure  economical,  the  choice  is  very  natural,  since  pulses 
are  instrumental  in  the  definition  of  integration  [16,  Chapter  10], 

This  choice  has  proven  very  successful,  nevertheless,  as  is  evident 
from  the  performance  of  the  solution. 

Circular  posts  and  thin  diaphragms  and  windows  were  treated 
early  using  the  Variational  Method  [1,  Chapter  3],  [2,  Section  8-9], 

In  [18,  Section  6-3],  the  "Singular  Integral  Equation  Method”  was 
applied  to  the  symmetrical  thin  window.  In  contrast  to  the  limited 
application  of  these  two  methods,  the  present  analysis  is  quite  general. 
It  can  become  the  first  step  in  the  solution  of  the  system  of  multiple 
dielectric  posts  in  the  capacitive  position  in  a  rectangular  waveguide. 
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Appendix  A 


Consider  the  function  defined  by  the  series 


Gi =  l 

n=l 


cos(nri)  cosCnn1)  e 

n 


-no  ^ 


o  -  k-c'l  • 


Since 


cos(nn)  cos(nr)1)  =  ~x  (cos  n(r)-r|!)  +  cos  n(r]+r)f)) 


(A.l) 


\  Re  (e-jn(rl-n,)  +  e-J-'l+n')) 


(A. 2) 


(A.l)  becomes 


G1  "  2  Re(  \  n 
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1  e-n(a+j (n-n')> 


+  l  -e 

n=l  n 


l  -n(o+j (n+n')) 


)  •  (A. 3) 


In  (A. 2)  and  (A. 3),  Re(z)  denotes  the  real  part  of  z, 
Since 


uo 

1  v  n 

T=i  =  +  2  * 

n=0 


z  e  C,  z  <1 


(A.  4) 


and  the  series  in  (A. 4)  converges  uniformly  for  all  z,  |z|  5  lz0l 
a  term  by  term  integration  can  be  carried  out  [20,  Section  5-4],  giv¬ 
ing 
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(A. 5) 
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In  (A. 5),  log  denotes  the  natural  logarithm.  Putting  equal  to 
e  ^  ^  ^  and  e  (n+H  ))  (^,5),  then  using  the  results 

in  (A. 3) ,  is  found  to  be 

Gl  -  -  i  Re [ iog((i  -  e-(°+3(n-n')))(1  .  e-<a+j(n+n,))))1 , 

Finally,  since 

|(1  _  e-(a+j(n-n,)))(1  _  e-(a+j(n+n,))) |2 

=  e  -  2e°  cos  (r]+n 1  )+l)  (e^°  -  2e°  cos (rrr| !  )+l) 

=  4e  (cosh  a  -  cosCn+n*)) (cosh  a  -  cosCn-n1)) 

becomes 

Gn  =  4a-  T  log[4(cosh  a  -  cos(n+n'))  (cosh  a  -  cos(n-Ti'))]  . 
12  4 

Since 

a  -  c-i(^'>)(i  -  e-j(n-''')) 

=  i  _  e-j(n+n’)  _  e-j(n-n')  +  e~i 2n 

=  2e  '*r|  (cos  n  -  cos  n') 

G^,  as  (J  tends  to  zero,  becomes 

G.  =  -  ^  log  (2 | cos  n -  cos  n ' I )  • 

The  singular  part  of  can  readily  be  extracted  from  (A. 8). 


(A. 6) 


(A. 7) 


(A.  8) 


(A. 9) 


(A. 10) 


That  is,  clearly. 
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Gis  =  -  |  l°g (2  (cosh  (£-£')  -  cos  (n-rT))) 


(A. 11) 
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Chapter  5 
DISCUSSION 

Three  systems  of  waveguide  discontinuities  have  been  considered 
in  this  dissertation. 

The  first  system  is  that  of  multiple  apertures  of  arbitrary 
shape  in  the  transverse  plane  between  two  cylindrical  waveguides. 

The  second  system  consists  of  metallic  obstacles  in  a  rectangular 
waveguide  that  are  uniform  along  the  narrow  side  of  the  waveguide,  but 
are  otherwise  of  arbitrary  shape  and  thickness,  i.e.,  a  system  of  in¬ 
ductive  posts.  Finally,  the  third  system  consists  of  metallic  ob¬ 
stacles  in  a  rectangular  waveguide  that  are  uniform  along  the  broad 
side  of  the  waveguide,  but  are  otherwise  of  arbitrary  shape  and  thick¬ 
ness,  i.e.,  a  system  of  capacitive  posts.  Common  between  the  first 
and  second  systems  are  the  inductive  windows  and  strips  in  a  rec¬ 
tangular  waveguide,  and  between  the  first  and  third  systems  are  the 
capacitive  windows  and  strips  in  a  rectangular  waveguide. 

The  analysis  of  waveguide  discontinuities  has  primarily  been 
carried  out  on  an  elementary  scale.  That  is,  solution  techniques 
have  been  sought  and  applied  for  problems  with  a  single  disconti¬ 
nuity.  Rarely  have  problems  involving  three  or  more  discontinuities 
been  considered,  and  only  then  if  they  are  all  of  the  same  shape. 
Although  some  of  the  elementary  problems  are  of  practical  importance 
that  warrants  considering  them,  these  and  others  are  better  worked 
out  as  special  cases  of  general  systems  of  discontinuities.  This » 
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of  course,  provided  that  general  and  efficient  solutions  can  be  found 
for  such  systems.  It  has  been  the  purpose  of  this  dissertation  to 
demonstrate  that  for  the  three  systems  considered. 

Other  systems  of  discontinuities  can  be  handled  using  the  methods 
employed  in  this  dissertation.  In  fact,  the  solutions  given  in  Chap¬ 
ters  3  and  4  can  become  the  first  step  in  the  analysis  of  systems  of 
multiple  dielectric  posts  in  the  inductive  and  capacitive  position 
in  a  rectangular  waveguide,  respectively.  Another  important  system 
of  discontinuities,  but  one  for  which  major  revisions  in  the  analysis 
of  Chapter  3  are  to  be  made,  is  that  of  metallic  resonant  posts  in 
the  inductive  position  in  a  rectangular  waveguide.  These  systems,  as 
well  as  others  that  can  be  identified,  are  recomended  for  future 
study . 

It  is  believed  that  advances  can  be  made  in  the  study  of  the 
waveguide  discontinuities  by  considering  systems  of  discontinuities 
rather  than  individual  discontinuities.  General  and  efficients  solu¬ 
tions  can  be  developed  for  these  systems  as,  it  is  hoped,  has  been 
demonstrated  in  this  dissertation. 
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